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摘 要 


高 速 侵 彻 数值 模拟 研究 通常 采用 拉 格 朗 日 有 限 元 法 。 < 
We- 适当 的 
人 
流体 动力 学 (SPH) 方法 模拟 5.56 mm M193 步枪 弹 侵 彻 钢板 过 程 ， 研 究 单 元 类 型 与 仿 
真 算法 对 模拟 精度 的 影响 ,分析 铬 芯 弹 穿 甲 过 程 和 原理 ,并 与 实验 结果 比较 。 结 果 表 明 ， 
SPH 法 在 模拟 铅 芯 弹 高 速 冲 击 问 题 中 具有 最 高 的 仿真 精度 、 稳 定性 和 效率 ， 可 准确 复 现 
实验 中 靶 板 的 变形 、 侵 蚀 、 剖 塞 过 程 。 高 速 铅 心 弹 对 钢 划 的 破坏 主要 分 为 半 流 体 侵 彻 、 
背 表面 损伤 、 剪 切 冲 塞 三 部 分 : 400 HBW 及 以 下 的 中 低 硬 度 装甲 钢 以 前 两 者 为 主 ，400 
HBW 部 由 于 拉 伸 塑性 较 低 而 易 发 生 背 面 月 落 ; 500 HBW 及 以 上 的 高 硬度 装甲 钢 以 后 者 
为 主 ， 界 面 击溃 发 生 于 $00~600 HBW 间 。 提 出 基于 Rosenberg 半 流 体 侵 彻 理论 和 Recht 
本 雪人 计 理 和 的 多 下 \ 弹 穿 甲 模型 ， 受 限于 模型 本 身 的 不 准确 性 和 弹 体 几何 的 复杂 性 ， 仅 
能 定性 描述 穿 甲 原理 ， 以 及 初步 预测 铅 芯 弹 对 靶 板 的 破坏 模式 。 


关键 词 : 光滑 粒子 流体 动力 学 ;有 限 元 仿真 ;高速 冲击 ; 半 流 体 侵 彻 
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Numerical Simulation of Lead Core Projectile [Impacting Steel Target 
at High Velocity 


Abstract 


Numerical simulation of high-velocity penetration ls usually carried out by Lagrangian 
finite element method. The process of lead-core 名 ll] metal jacket bullet Impacting steel target 
belongs to a typical large deformation and high strain rate problem，and the traditional 
Lagranglian method has poor Simulation effect on it. In order to find a suitable Simulation 
method, this paper adopts Lagranglian finite element method (FEM)，element-free Galerkin 
(EFG) and smoothed particle hydrodynamics (SPH) to simulate the process of 9.36 mm M193 
penetrating steel target studies the Influence of element form and simulation algorithm on the 
accuracy,，analyses the penetration process and principle of lead-core bullet， and compares 
with the experimental results. Results show that SPH has the best accuracy，gstability and 
efficiency in _ Simulating the high-velocity impact problem of lead-core bullet，and can 
accurately reproduce the deformation，eroslion and plugging process of the target plate in 
experiment. The damage of steel target by high-velocity lead-core bullet can be divided into 
three parts: Semi-fluid penetration, back-surface chipping and shear plugging. The former two 
are dominant for medium and low hardness armour steel of 400 HBW and below, and the 
back-surface chipping ls easy to occur in the target of 400 HBW due to its low tensile 
plasticity; whpile the latter ls dominant for high hardness armour steel. A penetration model of 
lead-core bullet based on Rosenberg Semi-fluid penetration theory and Recht shear plugging 
theory is proposed. Limited by the inaccuracy of the model itself and the complexity of the 
projectile geometry, it can only describe the penetration principle qualitatively and predict the 
damage pattern of lead-core bullets on the target plate preliminarily. 


了 Key Words: Smoothed Particle Hydrodynamics; Finite 上 Element Method; High-velocity 


Jmpact; Semi-fluid penetration 
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引 


了 路 


铅 芯 一 铀 被 弹 〈Lead-core Full Metal Jacket) 是 最 典型 的 一 种 步枪 弹 结 构 ， 其 起 源 可 
追溯 至 19 世纪 末期 的 瑞士 的 施 密 特 一 鲁 宾 步枪 ， 旨 在 提高 枪弹 整体 强度 以 减少 在 内 弹 
道 过 程 中 的 弹 体 变 形 和 枪 管 磨损 。 随 着 单 兵 及 轻型 车 辆 防护 技术 的 发 展 ， FMJ 穿 甲 威力 
弱 的 劣势 愈 发 明显 ， 故 正 逐渐 被 钢 忌 弹 奉 代 册 。 然 而 ， 由 于 具有 杀伤 力 强 赔 、 射 击 精 
度 高 B 的 优点 ， 时 至 今日 ，FMJ 结构 仍 被 广泛 应 用 于 下 至 小 口径 步枪 、 上 至 高 射 机 枪 的 
各 类 步 /机 枪弹 中 。 除 此 以 外 ， 铅 芯 弹 穿 甲 性 能 较 钢 一 铅 芯 组 合 的 半 穿 甲 弹 稳定 内， 因此 
也 被 大 量 用 于 评估 轻型 装甲 的 抗 弹 能 力 口 。 

要 对 目标 造成 有 效 破坏 ， 枪 弹 必 须 具 有 足够 的 威力 以 击 穿 相应 的 防护 装备 。 铅 蕊 弹 
的 主要 杀伤 对 象 为 装备 有 防弹 衣 等 的 有 生 目 标 , 其 防护 设备 材料 通常 为 钢 、 纤 维和 陶 盗 。 
早期 单 兵 防护 设备 通常 由 高 强度 钢板 冲压 而 成 ， 厚 度 一 般 在 2 mm 以 下 ， 仅 能 有 效 阻 挡 
手枪 弹 。 随 着 成 分 设计 、 治 炼 工艺 和 热处理 工艺 的 不 断 发 展 ， 高 性 能 装甲 钢 不 断 涌现 。 
从 中 低 强 度 向 高 硬度 、 超 高 硬度 、 高 强度 和 高 韧性 发 展 铅 。 当 代 装 甲 钢 依 和 硬度 分 为 均 质 
轧 制 装甲 《Rolled Homogeneous Armour) 、 高 硬度 装甲 钢 (High-Hardness Armour) 和 
超 高 硬度 装甲 钢 (Ultra-High-Hardness Armour) 三 类 。RHA 硬度 在 400HBW 以 下 ， 在 
保证 综合 抗 弹 能 力 的 同时 兼 具 良好 的 焊接 、 冷 弯 和 切割 性 能 ， 主 要 用 于 制造 几何 形状 复 
杂 且 厚度 较 大 的 坦克 装甲 四，HHA 硬度 约 为 S00 HBW， 兼 顾 硬度 和 韧性 ， 广 泛 用 于 制 
造 各 类 单 兵 防护 设备 〈 如 防弹 衣 ) 鲜 和 坦克 、 装 甲 车 等 的 复合 装甲 ，UHHA 硬度 在 600 
HBW 以 上 ， 多 用 于 轻型 装甲 车 的 披挂 装甲 。 为 综合 体现 铅 心 弹 对 各 类 目标 的 穿 透 能 

本 文 所 用 靶 板 硬度 在 300~600 HBW 间 ， 帮 盖 现代 装甲 钢 的 主要 硬度 范围 。 

目前 枪弹 穿 甲 研究 以 高 硬度 钢 芯 的 穿甲弹 为 主 ， 对 铅 芯 弹 的 探讨 相对 较 少 ， 且 多 局 
限于 实验 。Kile 等 站 通过 实验 与 数值 模拟 研究 了 7.62 mm 钢 芯 穿 甲 弹 对 高 强度 钢 靶 的 侵 
彻 过 程 ， 发 现 SPH 模拟 的 刘 板 对 弹 体 的 阻力 较 FEM 法 人 往 低 ， 但 能 更 好 地 表现 骨 落 和 碎 
片 飞溅 过 程 。 然 而 SPH 存在 多 种 算法 ， 以 LS-DYNA 为 例 ， 其 可 选 近似 方法 总 计 多 达 
16 种 ， 而 原文 仅 使 用 默认 方法 ， 未 做 进一步 探究 。 李 勇 9 通 过 数值 模拟 研究 14.5 mm 
钨 合 金 罕 甲 弹 对 钢 部 的 侵 彻 过程 ， 发 现 当 使 用 对 称 边 界 时 ，FEM、SPH、SPH-FEM 耦 
合 三 种 方法 计算 结果 相差 很 小 ， 均 能 准确 复 现 延 性 穿 甲 产 生 的 花 为 形 翻 唇 ; SPH 的 计算 
耗 时 高 于 FEM。Bervik 等 0 研究 了 7.62 mm 铅 芯 弹 和 钢 芯 弹 对 五 种 高 强度 钢 的 侵 彻 过 
程 ， 发 现 单 屋 和 无 间 辽 的 双 层 靶 具 有 相同 的 抗 弹 性 能 ， 而 有 间 隐 的 双 层 靶 劣 于 前 两 者 ; 

维 轴 对 称 的 FEM 和 多 物质 ALE 都 不 能 可 靠 地 模拟 铅 芯 弹 的 侵 彻 过 程 。Magier 等 上 2] 
使 用 单 物 质 ALE-FEM 耦合 法 模拟 7.92 mm DS 反 坦 克 步 枪弹 对 早期 装甲 钢 的 侵 彻 过 程 ， 
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与 历史 数据 对 照 具 有 良好 的 精度 ;并 发 现 高 速 下 铅 世 与 靶 板 会 相互 侵蚀 ， 铬 芯 与 被 甲 的 
碎片 沿 孔 壁 反 向 流动 ， 过 程 类 似 于 高 速 长 杆 弹 的 半 流 体 侵 彻 。Coghe 等 5 通过 实验 获得 
了 5.56 mm M193 铅 芯 步 检 弹 在 不 同 速度 下 对 6.5 mm Secure 500 的 毁伤 形 犁 ， 发 现 铅 世 
弹 在 高 速 下 对 高 硬度 装甲 钢 的 破坏 模式 以 剪 切 冲 塞 为 主 ， 表 面 仅 发 生 少 量 侵蚀 ， 且 穿 甲 
能 力 较 前 钢 芯 、 后 铅 芯 的 SS109 更 优秀 ，Sherrill 等 0 通过 实验 发 现 较 低 速度 下 两 种 枪 
弹 表 现 与 高 速 下 相反 ， 且 铅 芯 弹 对 陶瓷 一 纤维 复合 装甲 穿 甲 效 能 显著 劣 于 半 钢 芯 弹 。 
Binar 等 03 通 过 实验 获得 了 M193 对 6 mm Armox 600 超 高 硬度 装甲 钢 的 弹道 极限 速度 。 
Siriphala 等 09 通 过 实验 获得 了 Ml193 对 7mm Armox 600 的 毁伤 形 瑶 , 并 发 现 FEM 无 法 
准确 模拟 侵 彻 过 程 。 可 见 ， 对 铅 心 弹 高 速 穿 甲 的 研究 仍 处 于 初步 阶段 ， 不 同 速度 、 不 同 
靶 板 强度 仿真 与 实验 所 得 结论 有 很 大 差别 ， 暂 不 存在 统一 、 广 泛 适 用 的 理论 模型 。 

在 基于 拉 格 衣 日 坐标 的 有 限 元 方法 〈finite element method) 中 ， 计 算 网 格 的 每 个 节 
点 在 其 运动 过 程 中 跟随 指定 的 粒子 ， 因 此 外 表面 和 接触 面 很 容易 被 追踪 。 由 于 标准 有 限 
元 方法 的 接触 算法 能 可 靠 定 义 ， 且 计算 较为 快速 ， 因 此 是 研究 侵 彻 问题 的 主要 数值 模拟 
工具 。 尽 管 具 有 以 上 优点 ， 拉 格 朗 日 有 限 元 法 对 大 变形 、 大 失效 工 况 一 一 如 本 文 研究 的 
高 速 铅 芯 弹 冲击 问题 一 一 处 理 不 佳 Q0023。 若 网 格 发 生 严重 扭曲 、 变 形 ， 则 可 能 导致 剪 切 
锁定 073; 或 由 于 单元 自身 节点 与 面相 互 穿 透 而 引发 负 体积 错误 ， 进 而 使 单元 在 发 生物 理 
失效 前 先 数值 失效 ， 以 致 抗 弹性 能 异常 偏 低 。 此 外 ， 由 于 网 格 法 必须 删除 失效 单元 ， 整 
个 计算 过 程 中 质量 和 动量 不 守恒 ， 故 不 适用 于 对 弹 计 惯 性 敏感 的 高 速 冲 击 问题 53; 且 由 
于 接触 界面 的 不 断 变化 而 引发 非 真实 的 接触 力 振荡 021， 进 而 可 能 导致 弹 靶 过 度 侵蚀 。 此 
外 ， 对 于 高 速 侵 彻 问题 ，FEM 还 存在 部 件 间接 触 穿 透 Pa 和 接触 力 偏 低 避 的 问题 。 无 网 
格 伽 辽 金 方法 (Element-Free Gelarkin) 部 分 摆脱 了 对 网 格 的 依赖 。 由 于 节点 间 不 存在 连 
接 ， 避 免 了 单元 扭曲 问题 ， 较 FEM 显著 提升 了 计算 稳定 性 。 然 而 尽管 EFG 被 归 类 为 无 
网 格 方法 , 但 依赖 背景 网 格 以 进行 数值 积分 , 因此 仍 存在 由 网 格 /单元 删除 引发 的 质量 不 
守恒 等 问题 。 与 其 它 无 网 格 方法 类 似 ，EFG 需要 对 边界 条 件 进行 特殊 处 理 24， 其 中 由 于 
与 标准 拉 格 朗 日 有 限 元 法 的 相似 性 ，FEM-EFG 耦合 被 认为 是 一 种 有 效 的 方法 。 

除 拉 格 朗 日 坐标 外 ， 欧 拉 坐 标 类 方法 由 于 对 大 变形 的 优秀 适应 性 也 在 数值 仿真 中 受 
到 广泛 应 用 。 目 前 在 高 速 侵 彻 问题 中 主要 应 用 的 是 任意 拉 格 朗 日 一 欧 拉 方 法 〈Arbitrary 
Lagrange-Euler) 。 但 ALE 法 的 材料 边界 不 清晰 ， 对 强度 模拟 也 不 准确 请， 因此 仅 适 用 
于 强度 较 小 且 大 变形 的 枪弹 ， 而 不 可 用 于 模拟 高 强度 、 小 变形 的 半 板 ， 故 必须 使 用 流 固 
耦合 (Fluid-Solid Interaction) 以 模拟 两 类 材料 的 交互 。 然 而 ， 所 有 ALE 方法 一 一 包括 
单 物 质 ALE、 多 物质 ALE 和 结构 化 ALE 一 一 在 与 固体 单元 耦合 时 均 存 在 流体 泄露 问题 ， 
必须 依 经 验 设置 参 数 才 可 确保 稳定 上 5， 无 法 被 用 于 工 况 复杂 多 变 的 大 规模 分 析 。 此 外 ， 
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由 于 ALE 是 Euler 方法 的 一 种 变 体 ， 其 单元 和 节点 在 整个 模拟 过 程 中 固定 、 材 料 在 共 节 
点 的 单元 间 流 动 ， 而 非 如 拉 格 朗 日 类 方法 一 般 网 格 和 材料 协同 运动 ， 无 法 直接 由 节点 准 
确 描述 部 件 的 形状 。 
光滑 粒子 流体 动力 学 法 〈Smooth Particle Hydrodynamics) 是 一 种 拉 格 明日 无 网 格 方 
法 。 通 过 一 组 插值 粒子 离散 系统 体积 ， 每 个 粒子 携带 着 计算 场 和 内 在 数据 的 离散 值 ， 并 
在 控制 方程 的 规律 下 移动 。SPH 作为 一 种 完全 的 无 网 格 方法 ， 能 从 根本 上 避免 网 格 扭曲 
带 来 的 一 系列 问题 ， 同 时 ， 由 于 采用 Lagrange 坐标 ， 也 不 存在 Euler 类 方法 难 定 位 的 缺 
点 。 因 此 ， 从 以 上 来 看 ，SPH 是 一 种 适用 于 模拟 铅 忌 弹 冲 击 钢 靶 的 方法 。 然 而 ， 更 广泛 
的 研究 表明 ，SPH 同样 存在 根本 性 缺陷 , 在 高 速 冲 击 中 主要 表现 为 应 力 不 稳 定 和 粒子 异 
第 穿 透 。 应 力 不 稳 定性 的 根本 原因 在 于 应 力 状态 与 核 函 数 的 不 匹配 : 负 压 状态 下 粒子 间 
产生 随 着 粒子 间距 的 减 小 而 增 大 的 引力 ， 导 致 拉 伸 不 稳定 性 ; 正 压 状 态 下 粒子 间 产 生 随 
着 粒子 间距 的 减 小 而 先 增 大 后 减 小 的 斥 力 ， 导 致 压缩 不 稳定 性 。 杨 秀峰 等 上 5 改进 了 传统 
SPH 方法 的 核 函 数 和 离散 格式 ， 从 原理 上 消除 了 此 问题 , 并 具有 理论 上 的 普 适 性 。 然 而 ， 
其 改进 方案 目前 尚未 经 过 工程 验证 ， 暂 无 法 确定 是 否 存在 其 它 方面 的 缺陷 。Xnu 等 eq 在 
LS-DYNA 的 基础 上 开发 了 上 自 适 应 、 重 归 一 化 和 Lagrange 近似 方法 ， 在 配合 四 阶 样 条 核 
函数 时 ， 三 者 均 能 较 好 地 抑制 了 拉 伸 不 稳定 的 问题 。Lagrange 近似 法 仅 可 被 用 于 求解 小 
变形 问题 ， 自 适应 核 函 数 需 搭配 使 用 椭圆 支撑 域 ， 其 计算 稳定 性 不 及 传统 球形 文 撑 域 ， 
在 发 生 大 量 失 效 时 此 缺陷 尤为 明显 ; 重 归 一 化 近似 广泛 适用 于 绝 大 多 数 固体 问题 ， 但 对 
拉 伸 不 稳定 的 抑制 效果 不 及 前 两 种 方法 。Liu 等 上 使 用 重 归 一 化 近似 方法 模拟 了 有 限 厚 
和 半 无 限 混凝土 靶 受 尖 头 弹 侵 彻 的 过 程 ， 准 确 复 现 了 实验 中 靶 板 正面 变形 与 靶 后 碎片 云 
的 形 貌 。Bojanowski 等 CE 研究 长 杆 弹 侵 彻 土壤 的 过 程 ， 发 现 SPH 和 多 物质 ALE 均 能 获 
得 准确 的 模仿 真 结 果 ， 其 中 SPH 能 与 FEM 良好 耦合 和 接触 ， 并 且 对 网 格 密度 和 时 间 步 
长 不 敏感 。 然 而 ， 以 上 对 SPH 的 应 用 研究 多 局 限于 刚体 侵 彻 ， 弹 体 使 用 无 变形 和 失效 
的 FEM 网 格 ， 其 结论 未 必 完 全 适用 于 接触 面 两 侧 均 存 在 大 量 失效 的 高 速 半 流 体 侵 彻 。 
综 上 所 述 ， 铬 芯 弹 高 速 冲击 /问题 在 理论 和 仿真 上 均 不 成 熟 。 本 文 利用 LS-DYNA 实 
现 FEM、EFG、SPH 三 种 方法 模拟 35.356X45 mm M193 步枪 弹 冲 击 钢 轩 的 过 程 ， 探 完 单 
元 算法 与 参数 设置 对 仿真 精度 的 影响 ,基于 以 上 所 得 模型 研究 铅 心 弹 侵 彻 不 同 强度 钢 驾 
板 的 过 程 ， 参 考 半 流体 侵 彻 理论 和 仿真 结果 提出 铅 芯 弹 穿 甲 公式 ， 为 揭示 铅 忌 弹 穿 甲 机 
里 和 针对 性 防护 设计 优化 提供 一 定 借 鉴 思 路 。 


并 
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1 数值 模型 建立 


1.1 有 限 元 模型 


由 于 SPH 方法 对 边界 处 理 不 佳 ， 且 计算 耗 时 远 高 于 传统 网 格 法 PCJ， 因 此 本 文中 所 有 
SPH 模型 实际 均 为 SPH-FEM 耦合 模型 , 仅 大 变形 .大 失效 的 弹 体 及 靶 板 中 心 区 采用 SPH 
粒子 单元 ; 变形 相对 较 小 ， 且 需 设 固定 边界 约束 的 靶 板 外 周 部 分 采用 六 面体 FEM 单元 。 
EFG 及 FEM 弹 体 采用 四 面体 以 准确 描述 其 较 复 杂 的 几何 外 形 和 避免 负 体积 错误 ， 部 板 
整体 采用 六 面体 单元 以 提高 计算 精度 。 所 有 模型 靶 板 中 心 区 及 弹 体 原始 网 格 边 长 均 为 
0.25 mm， 靶 板 外 周 网 格 由 内 向 外 逐渐 稀 疏 ， 最 大 边 长 不 超过 5 mm。SPH 部 件 间 使 用 粒 
子 近似 方法 自动 识别 接触 CCONT=0) ， 无 需 额 外 定义 ; SPH 靶 板 核心 区 SPH 粒子 与 外 
周 区 六 面体 单元 设置 固 连 点 面 接触 (CONTACT TIED NODES TO SURFACE) ， 其 中 中 
心 区 为 从 节点 ， 外 周 区 为 主 面 。FEM、EFG 将 弹 、 部 全 体 部 件 设 为 一 部 件 集 (SET PART 
LIST) ， 定 义 侵 蚀 单 面 接触 (CONTACT ERODING SINGLE SURFACE) ， 此 接触 方法 
会 自动 识别 所 有 接触 面 ， 避 免 碎 片 与 母体 相互 穿 透 。 此 外 ， 设 置 内 接触 CONTACT 
INTERIOR) 以 抑制 网 格 内 穿 透 。 接 触 搜 索 深 度 DEPTH 设置 为 $， 面 一 面 与 边 一 边 穿 透 
在 计算 时 均 被 检查 。 其 余 参数 保持 默认 。 

SPH 粒子 近似 方法 采用 重 整 化 近似 (FORM=1) ， 此 方法 适用 于 大 多 数 固体 问题 ; 
核 函 数 设 为 四 阶 样 条 曲线 (SPHKERN=1) ， 光 涓 长 度 (CSLH) 初 设 为 13， 此 组 合 能 
较 好 地 抑制 拉 伸 不 稳定 问题 ， 避 免 SPH 靶 板 核心 区 在 与 FEM 外 周 区 固 连 的 边界 发 生 剪 
切 数值 失效 。FEM 六 面体 单元 均 选 择 八 节点 常 应 力 单元 (ELFORM=1) ， 此 单元 较 全 
积分 单元 (ELFORM=2、-1、-2) 能 更 好 地 避免 因 负 体积 而 导致 的 数值 失效 ， 但 因 是 缩 
减 积分 而 存在 沙漏 问题 ， 故 在 沙漏 控制 (CONTROL HOURGLASS) 中 设置 体积 精确 的 
Flanagan-Belytschko 刚性 沙漏 控制 〈IHQ=$) 。EFG 单元 均 选 择 标 准 无 网 格 实体 函数 

CELFORM=41) ， 核 函数 的 归 一 化 扩张 参数 (DX、DY、DZ) 均 设 为 1.05， 边 界 处 理 
方法 在 两 组 模型 中 分 别 选 择 FEM-EFG 耦合 法 〈IEBT=3) 与 快速 转换 法 (IEBT=4) ， 
其 中 前 者 在 边界 处 理 上 较 可 靠 ， 后 者 具有 最 优 的 计算 效率 和 稳定 性 ， 分 别 记 为 EFG3、 
EFG4;， 领域 积分 方法 选择 改进 的 高 斯 积分 〈IDIM=3) 。 


(a) SPH 弹 体 (b)FEMVEFG 弹 体 
图 1.1 弹 体 模型 
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(al) SPH-FEM 辑 (b) 全 FEMFEM-EFG 加 
图 1.2 靶 板 模型 


1.2 材料 模型 

弹头 及 靶 板 材料 均 Johnson-Cook 本 构 模 型 和 Gruneisen 状态 方程 。Johnson-Cook 模 
型 广泛 应 用 于 机 械 加 工 、 爆 炸 与 高 速 冲 击 领域 ， 综 合 考 虑 了 材料 的 应 变 强 化 效应 和 应 变 
率 与 温度 对 强度 与 塑性 的 影响 。 其 中 应 力 一 应 变 关系 的 展开 形式 为 : 


洲 


忆 了 一 了 
CO =(4+Bs)+Chn 一 )[]- (一 一 一 )” 1.1 
CC4+5500+Cm 二 (六 六 (1) 


式 中 : 4、B、m、C、 六 均 为 材料 参数 ， 由 实验 确定 。4 为 准 静态 下 屈服 应 力 ; 下 为 应 变 
硬化 常数 ，7 为 应 变 硬化 指数 ; C 为 应 变 率 硬化 系数 ; mm 为 温度 软化 指数 。cer、se 分 别 
为 等 效应 力 和 应 变 ; 、 分 别 为 实时 和 参考 应 变 率 ; 从 思 . 公 分 别 为 实时 、 参 考 和 熔 
化 温度 ， 歼 取 室 温 22 ， 各 种 钢 熔点 统一 取 1537.9 'C。 

失效 应 变 sy 由 式 (1.2) 确 定 : 
E 了 -了 


9 之 (1.2) 


712 于 


2r=(LD+D exzp Dio )(L+D,n 


式 中 : Di~ D5 为 材料 参数 ， 由 实验 确定 。 o= 为 应 力 三 轴 度 ， 其 中 ar 为 静水 压力 。 


当 损伤 变量 刀 达 到 1 时 单元 发 生 失效 ， 其 计算 方法 为 ; 


志 。 了 Ac。 


(1.3) 


己 
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表 1.1 Johnson-Cook 强度 模型 主要 参数 
材料 4/MPa BMPa 27 C 11 
CuZn30P81 111.7 504.7 0.42 0.009 1.68 
铬 四] 1 63 0.163 0.41 1 
4340B0161] 719 456 0.093 0.008 1.03 
Mars 190B2 1193 500 0.676 0.00435 1.17 
Armox 500B3] 1470 702 0.199 0.00549 0.811 
Armox 600B4 1580 958 0.175 0.00877 0.712 
表 1.2 Johnson-Cook 失效 模型 主要 参数 
材料 站 1 D2 1D3 2DD4 .1D5 
CuZn30BP3] 0 2.65 -0.62 0.28 0 
铅 P] 0 0 0 0 0 
4340B060 。 -0.8 2.1 -0.5 0.002 0.61 
Mars 190B2 0.21 7.2 -5.44 0 0 
Armox 500B3] 0.068 5.382 -2.554 0 0 
Armox 60064 -0.4 1.5 -0.5 0 0 


1.3 仿真 可 靠 性 验证 


实验 表明 , 6 mm Armox 600 钢板 可 保 误 


阻挡 937 mys 的 M19365], 弹道 极限 速度 (J 六 o) 


约 为 990 m/sn5。 由 于 Armox 600 钢板 存在 0.0/+0.6 mm 制造 公差 63， 故 可 认为 当 子 弹 
速度 为 990 m/s、 靶 板 厚 为 6.0 mm 时 能 恰好 击 穿 ， 剩 余 速 度 、 动 能 接近 于 0。 采 用 FEM、 
EFG3、EFG4、SPH 四 种 方法 模拟 以 上 两 种 速度 工 况 ， 残 余 质 量 (mr) 、 速 度 〈 万 ) 、 

动能 〈Eur) 、 开 坑 深 度 关 与 最 大 直径 Daax 如 表 1.3 所 示 。 其 中 SPH 法 穿 出 体 由 靶 板 与 


弹 体 的 碎片 组 成 ， 两 部 分 粒子 的 速度 差异 很 大 。 昌 


于 穿 出 体 以 塞 块 为 主 ， 故 速度 由 其 顶 


部 上 随机 4 节点 轴 疝 速度 取 平 均值 近似 求 得 。 因 仅 有 少量 弹 体 穿 出 ， 而 反 向 、 径 向 飞散 


的 弹 体 碎 片 也 共有 


定 动 能 ， 知 统计 在 内 会 产生 很 大 误差 。 因 此 残余 动能 以 丢 体 动能 近 


似 估计 。 考 虑 到 驾 体 的 整体 运动 ，Eur 实际 值 应 当 低 于 近似 结果 。 残 余 质量 以 残余 动能 
与 速度 估算 。SPH 方法 所 得 靶 板 抗 弹 能 力 最 低 ， 击 穿 、 阻 挡 情 况 与 实验 相同 053， 毁 伤 形 
钥 与 实验 结果 最 接近 09: 靶 板 受训 塞 穿孔 破坏 ， 冲 击 面 仅 受 轻微 侵蚀 。FEM 法 所 得 击 


穿 、 阻 挡 情 况 与 实验 结果 相同 ， 毁 伤 形 狐 与 实验 引 
锥 形 凹 坑 ， 凹 坑 底 部 发 生 延 性 罕 孔 。 两 种 EFG 法 所 得 驾 板 抗 弹 能 
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才 果 相差 很 大 : 靶 板 受 明 显 侵蚀 ， 形 成 


最 高 ， 击 穿 、 阻 挡 
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情况 与 实验 结果 不 一 致 ， 毁 伤 形 貌 与 实验 结果 相差 较 大 大 : 靶 板 受 明 显 侵蚀 ，EFG3 形 
成 半 柱 形 一 半 锥 形 凹 坑 ，EFG4 形成 柱 形 凹 坑 。 


表 1.3 四 种 仿真 方法 比较 
仿真 方法 vomsl mg wmsl 可 Doymm Danavymm 
FEM 937 ”未 击 穿 ， 靶 板 背 面 无 损伤 7.1 7.5 
FEM 990 ”0.04 241 2.51 ， 击 穿 7.3 
EFG3 937 “未 击 穿 ， 部 板 背 面 无 损伤 ”4.1 7.0 


恰 


EFG3 990 击 穿 ， 误 板 背面 开裂 6.1 6.8 
EFG4 937 ”未 击 穿 ， 部 板 背 面 无 损伤 2.8 7.9 
EFG4 990 “未 击 穿 ， 靶 板 背 面 无 损伤 ”4.1 7.8 
SPH 937 ”未 击 穿 ， 靶 板 背 面 无 损伤 ”3.1 11.3 
SPH 990 7.0 186 121  ， 击 穿 9.6 


(a) FEM (b) EFG3 (c) EFG4 (d) SPH 
图 1.3 990 m/s M193 侵 彻 6mm Armox 600 终 时 刻 弹 靶 Johnson-Cook 损伤 


对 比 四 种 仿真 模型 在 990 my/s 时 靶 板 毁伤 形 钥 与 实验 参考 ， 可 见 以 Armox 600 为 代 
表 的 超 高 硬度 装甲 钢 在 受 铅 忌 弹 冲 击 时 ， 靶 板 损 伤 分 为 挤 压 侵蚀 和 剪 切 冲 塞 两 部 分 ， 实 
验 和 SPH 以 前 切 训 塞 为 主 。SPH 粒子 在 达到 物理 失效 后 应 力 被 归 零 并 失去 体积 响应 ， 
但 仍 存在 于 积分 域 中 ， 且 保持 与 未 失效 粒子 的 接触 。 失 效 粒子 在 残余 弹 体 与 间 板 间 形 成 
了 一 层 无 应 力 、 几 乎 无 厚度 的 缓冲 层 ， 避 免 了 因 单 元 删除 而 产生 的 非 真实 空 除 ， 单 元 应 
力 不 会 因 接触 的 突然 断 开 和 恢复 而 振荡 。SPH 粒子 间 的 自动 接触 较 网 格 单元 间 的 接触 更 
平滑 ， 接 触 面 上 压力 更 均 义 ， 减 少 因 非 真实 集中 应 力 导 致 异常 失效 。 以 上 两 点 使 SPH 
刘 受 挤 压 侵 蚀 更 少 。FEM 及 EFG 法 计算 受 限于 网 格 ， 失 效 单元 在 仿真 过 程 中 被 删除 ， 
因而 不 存在 无 应 力 绥 冲 层 ， 刘 板 的 抗 挤 压 侵 蚀 能 力 劣 于 实际 情况 。 与 以 压缩 破坏 为 主 的 
长 杆 弹 侵 彻 土壤 2 类似，EFG 靶 抗 压 优 于 FEM， 其 中 EFG3 由 于 单元 行为 类 似 于 FEM 
其 挤 压 破 坏 程度 也 介 于 FEM 与 EFG4 之 间 。 


IE 
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从 图 1.4 阻力 变化 过 程 可 以 看 出 ， 与 前 文 对 靶 板 侵蚀 分 析 类 似 ，SPH 粒子 失效 前 后 
速度 和 动量 不 变 B9， 始 终 保持 对 靶 板 的 神 击 ,因此 铅 芯 和 整个 弹 体 所 受阻 力 均 明 显 高 于 
FEM 和 EFG。 在 冲击 面 直 径 与 其 它 方法 相差 不 大 的 前 提 下 ，SPH 部 的 环形 剪 切 面 分 配 
到 更 大 的 剪 应力。 由 于 侵 彻 过 程 中 弹 体 截面 积 和 肢 板 剩余 厚度 均 会 不 断 变 化 ， 无 法 直接 


200 了 


二 
加 


Resistance (kN) 


确定 其 每 个 时 刻 的 具体 值 ， 故 表 1.4 取 阻 力 峰 值 时 刻 近似 作为 最 大 前 应 力 时 刻 分 析 ， 其 
中 玉 为 最 大 阻力 ， 记 为 实时 接触 面 直径 ，! 为 靶 板 剩余 厚度 ，r 为 实时 至 


党 全 EFG4 Core 


和 生生 EEFG4Total 
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Time (hs) 


1.4 990 m/s M193 侵 彻 6mm Armox 600 过 程 阻 力 


表 1.4 990 my/s 侵 彻 仿真 中 的 冲 塞 剪 切面 


仿真 方法 FAN DAmmm bmm YVMPa 
FEM S$ 0.8 .2 4935 


EFG3 71 7.4 3.7 825 
EFG4 85 7.6 4.1 8068 


SPH 209 8.0 4.9 1379 


论 前 应 力 。 


根据 Johnson-Cook 横 型 , 在 常 光 准 静 态 条 件 下 , Armox 600 的 抗 剪 强度 为 1227 MPa， 


FEM、EFG3、EFG4 所 得 最 大 到 


最 大 到 


论 剪 应 力 均 显 著 低 于 此 值 ， 故 不 发 生 剪 切 破坏 ; SPH 的 


论 剪 应 力 高 于 此 值 ， 剪 切面 括 裂 并 形成 近似 柱 形 塞 块 。 理 论 分 析 与 图 1.3 中 的 仿 


真 结果 相 比 可 见 ，EFG3、EFG4 和 SPH 所 得 的 剪 切面 均 不 是 假设 的 规则 圆柱 面 ， 而 是 


受 背 表面 弯曲 产生 的 拉 应 力 与 正 表 面 训 击 产生 的 压 应 力 共 同 影响 、 直 径 先 增 大 后 减 小 的 
双 锥 组 合 面 ， 其 挤 压 端 直径 略 大 于 弹 坑 直径 。 但 从 图 1.6 (a) 中 的 Johnson-Cook 损伤 和 图 
1.3 (d) 的 最 终 冲 塞 孔 来 看 ， 综 合 应 力 与 应 力 三 轴 度 的 影响 ， 实 际 受 破坏 的 仍 是 与 理论 假 
设 尺 寸 、 形 状 均 相 近 的 圆柱 面 。FEM 法 挤 压 面 侵蚀 严重 ， 呈 近似 半球 形 ， 受 压 区 沿 径 
向 扩张 明显 ， 因 此 剪 切面 仅 有 锥 形 部 分 。 


_8- 


铅 芯 弹 高 速 冲击 钢 靶 仿真 研究 


Triaxialiy Factor (pivm) 
3.300e.01 


Triaxialiy Factor (pkvm) Triaxialiy Factor (pkvm) 


23 和 3.300e.01 
2.640e01 2.640e01 2.640e.01 
1.980e-01 1.980e-01 1980e.01 
1.320e-01 1.320e-01 1320e.01 
4 6.600e.02 
6.600e-02 e- 6.600e.02 
-2.776e-17 -2.776e-17 
-2.776e-17 
-6.600e-02 -6.600e-02 
-1320e01 -1320e.01 人 
-980e.01 -1.980e.01 -1320e41 
2.640e01 -2.640e.01 -1.980e.01 
-3.300e.01 -3.300e.01 .| -2.640e-01 
-3.300e.01 | 
人 
攻 学 
Ts 扣 寺 3 了 


图 1.$ 990 m/s M193 冲击 6 mm Armox 600 最 大 阻力 时 刻 靶 板 应 力 三 轴 度 


高 速 神 击 下 靶 板 受 应 变 率 强化 和 温度 软化 效应 影响 ， 实 际 抗 剪 应 力 随 之 增 减 。 以 与 
实验 结果 最 接近 的 SPH 法 为 例 〈 图 1.6 (bD)) ， 剪 切面 上 的 Tresca 应 力 约 为 1271 MPa， 
表明 应 变 率 强化 影响 略 大 ， 但 与 准 静态 一 常温 值 相 比 增 幅 很 小 ， 初 步 分 析 时 仍 可 用 后 者 
近似 估算 。 虽 然 理论 最 大 剪 应 力 可 用 于 粗略 判断 是 否 发 生 冲 塞 ， 但 由 于 剪 切 破坏 并 非 仅 
与 强度 有 关 的 瞬间 过 程 ， 因 此 分 析 准 确 性 较 差 。 此 外 ，FEM 和 EFG 虽 在 较 小 变形 情况 
下 计算 效率 较 高 CJ， 但 由 于 网 格 畸变 会 导致 的 时 间 步 缩短 ， 且 为 避免 单元 负 体积 错误 和 
网 格 异 常 穿 透 需 缩小 时 间 步 长 系数 (TSSFAC ) 以 确保 稳定 性 , 实际 计算 时 长 均 高 于 SPH。 
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(a]) ”Johnson-Cook 损伤 (b)  Tresca 最 大 剪 应力 
图 1.6 SPH 法 模拟 990 m/s M193 冲击 6mm Armox 600 15 ns 时 刻 刘 板 破坏 情况 
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1.4 本 章 小 结 

本 章 建立 了 FEM、EFG3、EFG4、SPH 铬 芯 弹 高 速 冲击 仿真 模型 ， 分 析 四 种 单元 算 
法 的 特点 。 结 果 表 明 ，SPH 法 具有 最 佳 的 计算 效率 、 精 度 和 稳定 性 ， 能 准确 复 现 靶 板 的 
剪 切 冲 塞 破 坏 ， 因 此 以 下 仿真 均 采 用 SPH 法 进行 。 


二 二 人 = 
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2 SPH 仿真 参数 优化 


2.1 SPH 方法 基本 原理 

SPH 方法 最 初 由 工 . B. LucyB1 和 开 . A. Gingold65 提 出 ， 用 于 求解 天 体 物 理 问题 。 此 
方法 的 进一步 发 展 使 得 它 的 实用 性 也 扩展 到 了 固体 和 流体 力学 的 广泛 问题 上 。 为 表征 
SPH 粒子 在 时 间 步 内 的 运动 ， 需 要 构造 一 个 近似 函数 ， 其 基本 形式 为 ; 


JD= 忆 十 /CD 太一 胃 C. 
式 中 下 、 户 为 某 一 附近 点 j 的 质量 和 密度 ， 球 为 核 函 数 ， 其 求解 方式 为 ， 
Ji 一 必 ) 
矿 (x， 9 (2.2) 


式 中 4 为 空间 维 数 ， 本 文 所 有 模型 均 为 三 维 ，20(Cz) 为 辅助 机 数 ， 如 五 阶 样 条 函数 ， 由 参 
数 SPHKERN 控制 .7 为 光滑 长 度 , 用 于 判断 附近 点 j 对 目标 点 站 的 影响 程度 。LS-DYNA 
中 的 SPH 求解 器 使 用 可 变 的 光滑 长 度 。 通 过 各 部 件 中 节点 间距 的 最 大 、 最 小 值 和 CSLH 
常数 计算 初始 光滑 长 度 zmo， 其 中 CSLH 推荐 值 为 1.2， 以 保证 计算 域 在 各 方向 上 均 能 履 
盖 两 层 以 上 粒子 。CSLH 越 大, 计算 域 覆 盖 面 积 越 大 。 若 值 小 于 1.05 会 导致 计算 不 稳定 ， 
大 于 1.3 会 显著 增加 计算 耗 时 。 各 粒子 具有 随时 间 变 化 的 、 独 立 的 光滑 长 度 ， 其 计算 方 
法 为 : 


d 1 
了 [ADI= 了 ADYY (2.3) 


式 中 V.v 为 流动 发 散 性 ;AD 为 光滑 长 度 ， 应 介 于 最 小 值 HMINX (Ci=0) 和 HMAXX 
CE=0) 之 间 。HMIN 默认 值 为 0.2。HMAX 与 核 函 数 有 关 : 对 于 三 阶 样 条 函数 ， 默 认 值 
为 2.0; 对 于 五 阶 样 条 函数 , 推荐 取 3.0 或 更 大 。 以 上 参数 中 , 可 在 LS-DYNA 中 选择 的 、 
对 仿真 结果 影响 最 明显 的 是 CSLH。 此 外 还 包括 粒子 近似 理论 FORM 和 失效 粒子 控制 方 
法 了 EROD。 以 下 通过 仿真 研究 此 三 参数 对 仿真 精度 的 影响 ， 进 而 确定 最 佳 参数 组 合 。 
2.2 初始 光滑 长 度 

首先 测试 初始 光滑 长 度 参 数 CSLH 的 影响 ， 网 格 /粒子 模型 及 材料 参数 与 前 文 相同 。 
分 别 对 初速 937 m/s 和 990 mys 的 侵 彻 过 程 进行 仿真 ， 以 击 穿 / 阻 挡 情 况 、 穿 孔 和 塞 块 形 
瑶 、 全 弹 体 和 铅 蕊 所 受阻 力 评 估 其 计算 准确 性 。 参 考 推荐 值 ，CSLH 分 别 设 为 1.05、1.1、 
1.2、1.3、1.4。 两 组 CSLH=1.05 均 发 生 粒 子 超速 错误 导致 计算 异常 终止 ， 以 下 不 再 分 析 。 
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由 图 2.1 可 见 ， 弹 体 所 受阻 力 随 CSLH 增加 而 增加 ， 均 呈 增 加 一 稳定 一 衰减 趋势 ， 主 要 
差距 在 于 稳定 段 的 持续 时 间 ，CSLH 越 低 ， 稳 定 阶段 越 短 。 由 于 SPH 法 存在 拉 伸 不 稳定 
问题 ,间距 过 大 的 粒子 会 在 发 生 由 Johnson-Cook 模型 控制 的 物理 失效 前 先 断 开 连 接 ,， 即 
数值 失效 。CSLH 1.1 靶 板 背面 受 拉 伸 破 坏 较 重 ， 划 板 在 被 冲 出 塞 块 前 先 数值 失效 产生 
小 缺口 ， 而 后 剩余 粒子 受 拉 应 力 破 坏 失 效 ， 与 实验 结果 有 明显 偏差 ; 在 阻力 上 表现 为 无 
稳定 侵 彻 阶段 ， 衰 减 段 阻力 下 降 更 早 。 此 外 ，CSLH 1.1 也 唯一 在 937 my/s 时 被 击 穿 : 与 
990 myAs 类 似 , 靶 板 背面 撕 裂 骨 落 产生 雁 片 云 ,弹头 未 穿 出 , 因此 显然 是 不 准确 的 。CSLH 
1.2~1.4 数值 失效 问题 逐渐 改善 ，1.3、1.4 阻力 曲线 几乎 重合 ;穿孔 均 为 近似 柱 形 ， 罕 出 
体 均 为 中 心 塞 块 + 环形 大 碎 块 + 弹 体 和 碎片， 与 6. mm Secure S00 FEM 仿真 和 实验 结果 接 
近 03]。 穿 孔 直 径 分 别 为 9.7 mm、9.6 mm、9.6 mm， 较 7mm Armox 600 实验 结果 略微 仿 
高 49。。 综 上 ， 可 以 认为 ，CSLH 1.3、1.4 均 较 好 地 抑制 了 拉 伸 不 稳定 及 其 引发 的 数值 
失效 问题 ， 能 准确 地 描述 阻力 过 程 。 虽 1.4 因 计算 时 间 较 长 而 不 在 推荐 范围 内 ， 但 对 上 毁 
伤 形 貌 和 阻力 过 程 描述 更 准确 ， 以 下 仿真 均 采 用 CSLH 1.4。 


二 中 二 CSLH 1.1 Core 
区 CSLH1.1 Total 
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(b) 990 mAs 
2.1 不 同 初始 光滑 长 度 模拟 不 同 速度 M193 侵 彻 6mm Armox 600 过 程 阻 力 
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三 直 业 二 


(a) CSLH=1.1 (b) CSLH=1.2 


2.3 粒子 近似 方法 
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(c) CSLH=1.3 
图 2.2 ”990 m/s M193 侵 彻 6mm Armox 600 不 同 CSLH 轩 板 毁伤 
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(0 CSLH=1.4 


LS-DYNA 内 置 了 多 种 粒子 近似 方法 , 其 中 可 用 于 SMP 计算 的 方法 名 称 与 编号 如 表 
4 所 示 。Scazzosi 等 B9 的 研究 表明 FORM 4、7、8 不 能 求解 高 速 冲击 问题 ， 故 在 此 不 作 
讨论 。 仿 真 结 果 显 示 ， 两 种 各 向 同 性 重 整 化 近似 方法 FORM 3、6 和 两 种 各 向 噶 性 近似 
方法 FORM 9、10 均 在 靶 板 SPH-FEM 边界 处 发 生 粒 子 超速 错误 ， 导 致 计算 无 法 继续 进 


行 ，SPH 与 FEM 无 法 通过 常规 的 固 连接 触 簿 合 ; 两 种 


曾 强 流体 方法 FORM 1$、16 均 有 


弹 体 碎片 粒子 发 生 超 速 错误 ， 导 致 计算 无 法 计算 进行 ， 并且 相 较 于 FORM 1 靶 板 背 表 面 
拉 伸 不 稳定 现象 更 严重 ; 所 有 非 重 整 化 近似 方法 FORM 0、2、5、15 均 与 图 2.3 (8) 类 似 ， 
在 靶 板 SPH-FEM 边界 处 发 生 数值 失效 ，SPH 与 FEM 部 分 分 离 ， 并 有 少量 边界 粒子 和 
实体 单元 发 生 异 常 物理 失效 。 此 外 ， 从 图 2.3 (b) 终 时 刻 靶 板 毁 伤 形 貌 来 看 ， 整 体 上 ， 相 


较 于 标准 方法 ， 重 整 化 近似 方法 能 更 好 地 抑制 拉 伸 不 稳定 问题 。 


一 致 ， 近 似 方法 选择 精度 最 高 的 重 整 化 近似 FORM 1。 


表 2.1LS-DYNA 


综 上 ， 以 下 与 前 文保 持 


可 用 于 高 速 冲击 问题 的 粒子 近似 方法 
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人 了 
(FORM0 (bj) FORMI1 
图 2.3 重 整 化 对 弹 、 间 形 貌 的 影响 


2.4 失效 粒子 控制 方法 

IEROD 关键 字 用 于 控制 失效 的 SPH 粒子 。IEROD 0 表示 粒子 始终 保持 活动 状态 ， 
其 应 力 不 会 在 失效 后 自动 归 零 ， 并 且 会 保留 体积 响应 ， 被 认为 不 适用 于 大 失效 问题 ; 
IEROD 1 表示 粒子 在 失效 后 部 分 失 活 ， 应 力 自 动 归 零 并 失去 体积 响应 ， 但 仍 保留 于 计算 
域 中 以 确保 计算 稳定 性 ;IEROD 2 表示 粒子 在 失效 后 完全 失 活 ， 应 力 自动 归 零 并 失去 体 
积 响 应 ;IEROD 3 保留 由 状态 方程 控制 的 体积 响应 ， 其 余 与 EROD 1] 相同。 仿真 发 现 ， 
IEROD 0 与 IEROD 3 在 FEM-SPH 边界 处 计算 稳定 性 很 差 , 在 使 用 多 种 粘性 和 刚性 沙漏 
室 制 〈IHQ=1、2、3、4、5) 时 均 因 在 边界 处 FEM 网 格 异 常 畸 变 发 生 负 体积 错误 进而 
导致 计算 终止 。IEROD 1] 与 EROD 2 以 半 流 体 侵 御 为 主 的 阻力 上 升 和 稳定 段 几 乎 重合 ， 
仅 下 降 段 后 期 有 明显 差别 ， 最 终 塞 块 动能 分 别 为 121、211 J。 与 IEROD 1 相 比 ，IEROD 
2 所 得 塞 块 主体 形状 基本 相同 , 但 以 拉 伸 破坏 为 主 的 头 部 边缘 相对 较 窄 ; 未 形成 碎片 环 ， 
而 是 一 个 较 窗 的 完整 环 与 大 量 失 效 粒 子 云 。 从 结果 来 看 ， 压 缩 破坏 过 程 对 失效 粒子 控制 
方法 不 敏感 ， 主 要 影响 在 于 拉 伸 和 剪 切 破坏 ，IEROD 2 靶 板 的 抗 拉 、 剪 韧性 相对 较 低 。 
此 外 ， 与 其 描述 相符 ，IEROD 2 计算 稳定 性 较 IEROD 1 略 差 ， 前 者 需 将 时 间 步 长 缩放 
系数 减 小 至 0.4 才能 稳定 计算 ， 否 则 会 发 生 粒 子 超速 错误 ， 而 后 者 仅 需 0.8。 综 上 ， 以 下 
仍 采 用 部 分 失 活 的 了 EROD 1 。 


2 IEROD1 


人 IEROD 2 
200 十 


3 
Time (hs) 


图 2.4 不 同 失 效 粒 子 控制 方法 模拟 990 m/s M193 侵 彻 6 mm Armox 600 阻力 过 程 
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History Variable 丰 
1.000e*+*00 
9.000e-01 
8.000e-01 
7.000e-01 
6.000e-01 
5.000e-01 
4.000e-01 
3.000e-01 
2.000e-01 
1.000e-01 
0.000e+*00 | 


Ca) IEROD 1 Cb) IEROD 2 
图 2.5 ”失效 粒子 控制 方法 对 靶 板 毁伤 形 貌 的 影响 


2.5 ”本章 小 结 


本 章 探究 了 初始 光 消 长 度 (CSLH) 、 粒 子 近 似 方法 (FORM) 和 失效 粒子 控制 方 
法 (IEROD) 对 仿真 的 影响 。 结 果 表 明 ，CSLH=1.4、FORM=1、IEROD=1! 的 参数 组 合 
能 最 准确 地 模拟 铅 芯 弹 冲 击 过 程 ， 且 有 具有 良好 的 计算 稳定 性 。 
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3 靶 板 硬度 对 穿 甲 过 程 的 影响 


3.1 靶 板 等 效 强度 


4340 的 性 能 与 其 热处理 方法 有 很 大 关系 , 不 同文 献 中 提出 的 数值 差别 较 大 。 由 于 文 

献 [30] 仅 给 出 了 洛 氏 硬度 (HRC) 、 必 服 强度 和 抗 拉 强度 ， 其 布 氏 人 硬度 依 GBMT 1172 换 
算 。 对 30 HRC0o 与 302 HBW00 钢 靶 抗 BS-41 侵 彻 的 研究 表明 ， 两 种 硬度 靶 板 的 抗 弹性 
能 不 存在 可 见 区 别 ,， 因 此 本 文 所 用 4340 可 近似 为 300 HBW 的 装甲 钢 。Mars 190 属于 典 
型 的 大 厚度 装甲 钢 ， 硬 度 随 板 厚 的 增加 而 显著 下 降 ， 本 文 及 文献 32] 均 取 薄 板 的 硬度 上 
限 。Armox 500 和 Armox 600 应 用 较为 广泛 ， 且 均 有 明确 的 生产 标准 ， 实 测 硬 度 分 别 为 
516 HBWI 和 5$79 HBW 上 AI]。 综 上 ， 尽 管 受 限于 生产 工艺 ， 各 部 的 硬度 都 存在 偏差 ， 但 
均 在 误差 允许 范围 内 。 
昌 于 装甲 钢 以 抗 弹 而 非 力 学 性 能 作为 指标 ， 其 强度 和 硬度 的 波动 范围 较 一 般 结构 钢 
偏 大 , 例 。 整体 来 看 , 屈服 强度 的 数值 约 相当 于 布 氏 硬度 CBHN) 的 3 倍 。 以 下 对 300~600 
HBW 范围 内 的 5 种 钢 靶 , 分 别 在 等 靶 厚 击 穿 、 临 界 击 穿 和 半 无 限 靶 三 种 工 况 下 进行 990 
m/s M193 侵 彻 仿真 ， 以 分 析 硬 度 对 钢板 抗 铅 蕊 弹性 能 的 影响 。 表 中 列 出 了 各 材料 的 名 
称 及 标 称 性 能 。Mars 190 较为 特殊 ， 其 便 弃 、 强 度 和 伸 长 率 范 围 极 大 ， 并 且 可 分 为 高 低 
两 级 硬度 的 等 级 ， 分 别 记 为 Class 1 和 Class 2。 与 其 硬度 、 功 能 相近 的 Armox 370 也 分 
为 两 级 ， 均 符合 MIL-A-12560 中 对 RHA 的 规定 。 本 文 使 用 的 是 硬度 较 高 的 一 双 


全 


表 3.1 靶 板 材料 标 称 性 角 
材料 PHBW omiVMPa omwMPa ， 伸 长 率 /% 
4340 300 850 1020 14~16 
Mars 190 Class 1 240~400 ”600~900 ”800~1000 12~20 
Armox 500 ， 480-540 1250 1450~1750 10 
Armox 600 570 - 640 ”1500 2000 7 


Johnson-Cook 模型 中 将 Taylor-Quinney 系数 设 为 0.9， 即 90% 的 塑性 功 转化 为 热 。 
应 变 率 在 高 速 侵 彻 中 影响 显著 ， 以 下 设 为 104/s。 图 3.1 (9) 展 示 了 等 效应 变 由 0 增长 至 2 
的 过 程 中 等 效应 变 的 变化 过 程 ， 可 见 四 种 钢 靶 的 变形 过 程 均 可 分 为 弹性 、 应 变 强化 和 热 
软化 三 个 阶段 ， 其 中 弹性 段 最 短 ， 应 变 强化 段 次 之 ， 热 软化 段 占 比 最 大 。 对 最 大 等 效应 
力 与 便 度 的 关系 进行 线性 拟 合 ， 结 果 如 式 (3.1) 所 示 ， 误 差 如 图 3.1 (b)。 
=3.528 玉 +80.54 (3.1) 


Oerr Imax 


二 | 本 二 
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忆 总 加 
CD Co 5 
CI 己 Cn 
] 所 二 


Stress (MPa) 
\ 
S 


(a) 应力 一 应 变 关系 


1 
分 
伺 
已 


Simulation 


Linear 


Maximum stress (MPa) 


Hardness (HBW) 
(b) “最 大 应 力 一 硬度 
图 3.1 考虑 应 变 率 强化 和 热 软 化 效应 的 应 力 一 应 变 关系 


钢 的 硬度 与 强度 的 关系 有 多 种 经 验算 法 。Tate 等 啊 认 为 在 长 杆 高 速 侵 彻 中 靶 板 材料 
可 近似 为 理想 弹 塑 性 ， 其 屈服 强度 orae 相当 于 4.2 倍 布 氏 硬 度 。JC 模型 线性 拟 合 的 截 
距 仅 相 当 于 orcmax 的 不 足 10%, 在 粗略 计算 时 可 近似 为 正比 。 但 Tate 式 的 斜率 明显 偏 高 ， 
以 致 过 度 高 估 了 装甲 钢 的 强度 。Rosenberg 等 以 准 静 态 一 常温 条 件 下 s0.27 时 的 应 力作 
为 等 效 强 度 沪 ]， 由 于 同样 久 视 了 热 软 化 ，oRes 在 数值 上 与 mmrate 接近 。 偏 差 可 能 来 自 于 
两 者 研究 对 象 的 不 同 : Tate 等 和 Rosenberg 等 所 用 邯 涵 盖 铝 、 铀 、 钢 等 各 种 金属 ， 硬 度 
普遍 在 300 HBW 以 下 ;而 本 文中 全 部 为 钢 靶 ， 且 孔 > 300HBW。 此 外 ， 各 金属 强度 
标准 均 存 在 硬度 一 抗 拉 强 度 插 值 换算 方法 ， 以 GBMT 1172 为 例 C9， 此 换算 方法 除 基 础 的 
硬度 外 还 考虑 了 钢材 化 学 成 分 的 影响 。 以 铬 镍 钢 估算 300 HBW 的 4340， 以 超 高 强度 钢 
估算 400~600 HBW 各 钢 靶 ， 结 果 如 om 所 示 。 尽 管 抗 拉 强 度 测 试 中 应 力 随 应 变 增 加 而 下 
降 的 原因 是 颈 缩 ， 而 与 热 软 化 无 关 ， 并 且 忽 略 了 应 变 率 强 化 效应 ， 但 两 者 呈现 出 良好 的 
相关 性 ， 各 钢 刘 的 偏差 均 在 10% 以 内 。 
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表 3.2 ”等 效 强度 估算 方法 比较 
等 效 强度 /MPa 误差 /% 
OJcmax OrTate aos am oft On 
300 1152 ”1260 1198 1046 1.1 9.2 
400 1450 ”1680 1448 1370 2.9 5.5 
500 1888 2100 2101 1716 2.3 9.1 
600 ， 2182 2520 2512 2180 0.7 0.! 


分 别 以 400 和 500 HBW 级 装甲 钢 作为 验证 :Hardox 400 的 orcmax 为 1484 MPar 7 
另 一 Armox 500 模型 为 1929 MPal4l, 30 PM 为 1993 MPat1, Protac 500 为 1991 MPab01。 
考虑 到 装甲 钢 的 硬度 和 强度 普遍 存在 一 定 生产 误差， 如 Armox 560T 标 称 最 小 屈服 强度 
1300 MPa51U， 实 测 1711 MPab55，Demox Protect 500 标 称 屈服 强度 1300 MPaB53， 实 测 
1592 MPab52; Hardox 400 标 称 屈服 强度 1100 MPa59， 实 测 1148 MPa52， 经 验 式 的 偏差 
在 合理 范围 内 。 综 上 ， 对 于 硬度 在 300~600 HBW 范围 内 的 装甲 钢 ， 可 由 Johnson-Cook 
模型 或 式 (3.0D) 估 算 其 在 侵 彻 过 程 中 的 最 大 应 力 ， 也 可 由 式 G3.D) 反 向 初步 评价 
Johnson-Cook 模型 的 可 靠 性 。 

为 研究 便 度 的 具体 影响 , 以 下 对 等 厚度 击 穿 、 弹 道 极 限 、 半 无 限 靶 三 种 工 况 下 M193 
的 侵 彻 进行 测试 ， 初 速 均 为 990 my/s。 


3.2 等 厚度 击 穿 


为 确保 各 组 均 能 被 击 穿 ， 取 和 靶 板 厚度 为 600 HBW 的 临界 击 穿 厚度 6 mm。 由 表 3.3 
可 见 ， 随 着 硬度 的 增加 ， 靶 板 的 破坏 形式 由 拉 伸 逐渐 过 渡 至 剪 切 ， 穿 透 后 弹 体 前 的 靶 板 
剩余 厚度 逐渐 增加 。300 和 400HBW 表现 为 拉 伸 破坏 ， 靶 后 雁 片 各 方向 分 布 较为 均匀 ， 
整体 上 附着 于 残余 弹 体 上 ， 呈 近似 半 权 球状 ; 500 和 600HBW 表现 为 剪 切 破坏 ， 其 中 
500 HBW 由 于 韧性 较 好 ， 剖 出 的 塞 块 较 600 HBW 时 完整 。 

依照 发 生 部 位 ， 靶 板 的 破坏 可 分 为 正面 压缩 、 环 形 前 切 和 背面 拉 伸 三 部 分 : 

正面 破坏 来 源 于 弹 对 部 的 直接 挤 压 ， 四 强度 靶 板 正面 破坏 均 仅 限 于 很 薄 的 一 层 ， 表 
明 在 6 mm 时 抗 压 不 是 抗 弹性 能 的 主要 因素 。 正 面 破坏 可 分 为 中 心 和 边缘 两 部 分 ;中心 
损伤 区 随便 度 的 提高 而 厚度 减 小 ; 边缘 损伤 区 随便 度 的 提高 而 增 大 。 尽 管 发 生 完全 破坏 
的 范围 很 小 ， 如 靶 体 损伤 一 栏 所 示 ， 但 类 似 于 长 杆 高 速 侵 彻 ， 塑 性 流动 也 会 造成 抗 阐 能 
力 的 损失 。 由 于 半 无 限 刘 高 度 承 压 ， 其 应 力 三 轴 度 很 低 ， 靶 板 的 失效 仅 局 限于 轴线 处 ， 


本 
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但 因 纯 粹 塑性 而 开 孔 。 因 此 , 以 侵 彻 结束 时 的 残余 靶 厚 评估 实际 正面 破坏 程度 , 其 中 300 
和 400 HBW 由 于 背面 破坏 较 强 

环形 破坏 类 似 于 典型 的 剪 切 种 塞 , 但 因 靶 板 存 在 变形 而 附加 了 拉 伸 的 影响 。 由 表 3.3 
中 靶 板 损伤 和 失效 两 行 可 见 ， 和 典型 的 环形 破坏 面 仅 存在 于 400~600 HBW。 随 着 靶 板 硬 
度 的 减 小 ， 剪 切 破坏 由 600 HBW 时 的 近似 柱 面 逐渐 扩张 为 400 HBW 的 约 45” 锥 面 ， 
同时 失效 面 的 厚度 逐渐 减 小 。300 HBW 由 于 受 正面 破坏 过 多 ， 坦 板 向 背面 显著 弯曲 ， 
因此 受 虽 没有 剪 切面 ， 但 由 临 穿 透 损伤 可 观察 到 拉 伸 破坏 。400 HBW 也 存在 相似 现象 。 

背面 破坏 来 源 于 凸 出 、 弯 曲 导 致 的 拉 伸 ， 分 为 表面 和 轴线 两 部 分 。 轴 线 破 坏 来 源 于 
弹 尖 的 集中 拉 伸 ， 表 面 破坏 来 源 于 整体 弯曲 。 刘 板 硬度 和 强度 越 高 ， 抗 弯 性 能 越 好 ， 同 
时 残余 靶 厚 越 大 ， 整 体 刚 度 越 高 背面 变形 越 小 ， 也 就 越 不 容易 发 生 拉 伸 破坏 。 也 与 高 应 
力 三 轴 度 下 的 塑性 有 关 。 例 如 对 于 400 HBW， 其 本 身 抗 弯 性 能 不 佳 ， 且 高 应 力 三 轴 度 
下 塑性 很 低 ， 因 此 背面 破坏 最 为 严重 ，300 HBW 虽然 抗 弯 强 度 不 及 400 HBW， 但 塑性 
更 好 ， 因 此 表面 破坏 较 轻 。 


表 3.3 弹 施 形 钥 
硬度 /HBW 300 400 500 600 


穿孔 形 狐 


临 穿 透 损伤 
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加 一 
加 
由 表 3.4 可 见 , 穿孔 直径 整体 上 与 靶 板 硬度 负 相 关 ， 主 要 原因 是 随 着 靶 板 硬度 提高 ， 
铅 亿 所 受 轴 疝 阻力 也 相应 增强 ， 进 而 铅 、 铜 径 向 流动 导致 的 扩张 越 强 。 


600 HBW 与 此 


趋势 异常 是 因为 靶 板 硬度 足够 高 ， 以 至 于 在 发 生 剪 切 冲 塞 前 几乎 不 受 侵蚀 ， 并 仅 有 轻微 


变形 。 


因此 最 终 的 穿孔 直径 更 接近 于 冲 塞 直径 ， 而 不 是 其 它 靶 的 侵蚀 直径 。 


残余 靶 厚 与 


辑 板 硬度 正 相 关 ， 即 靶 受 侵 es 与 硬度 负 相 关 。 其 中 300 HBW 虽然 弹 体 前 


方 没有 剩余 肢体 ， 但 部 并 非 被 完全 侵 


由 表 3.3 中 穿孔 形 租 一 栏 可 见 ， 因 为 残余 弹 体 
的 动能 较 高 且 头 部 较为 尖锐 ， ee 拉 断 。 两 侧 的 残余 靶 体 虽然 


RAY 


厚度 不 均匀 ， 但 从 图 像 可 定性 判断 薄 于 400 HBW， 即 与 前 文 所 述 趋势 相同 。 


表 3.4 990 m/s M193 临界 击 穿 数值 
便 度 /HBW 穿孔 直径 /mm ”残余 靶 厚 /mm 动 


300 8.88 
400 9.03 
500 9.53 
600 8.9 


3.3 弹道 极限 


由 表 3.4 可 见 ， 各 部 
的 环形 破坏 面 。 破 坏 面 由 


300 HBW 的 约 45”? 
表明 驾 板 在 冲 塞 前 发 生 的 弯 拉 变形 逐 


0 
1.03 
4.00 
4.60 


结果 

能 消耗 / 
1353 
1594 
1707 
1662 


三 2 人 


板 在 临界 击 穿 时 均 发 生 剪 切 冲 塞 ， 击 穿 损 伤 一 栏 中 均 可 见 明 下 
截 锥 形 逐 渐 收 缩 为 600 HBW 的 近似 柱 形 ， 
簿 渐 减 弱 。300 HBW 在 冲 塞 的 同时 发 生 与 6 mm 相似 


GE 
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的 拉 伸 破坏 ， 塞 块 与 靶 板 出 口 处 的 翻 展 被 拉 断 形成 大 碎 块 ; 400 HBW 各 在 高 应 力 三 划 
度 下 塑性 较 低 ， 发 生 背 面 朋 落 并 产生 完全 失效 的 碎 导 ， 以 致 塞 块 较 薄 ; 00 和 600 HBW 
为 典型 冲 塞 ， 其 中 S00 HBW 由 于 塑性 较 高 ， 塞 块 更 为 完整 。 各 间 临 界 击 穿 厚度 如 表 3.6 
中 7 Penmax 所 示 。 


表 3.5 临界 击 穿 工 况 下 弹 靶 形 狗 
300 HBVW 400 HBW S00 HBVW 600 HBVWV 


世 . 


弹 刘 形 貌 


驾 板 损伤 


击 穿 破 坏 


2 


包 
1 


> 
1 


习 
1 


全 


Maximum penetration thickness (mm) 


300 400 500 600 
Hardness (HBW) 


3.2 ”最 大 击 穿 靶 厚 


而 
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由 表 3.6 可 见 ， 背 面 变形 量 与 残余 半 厚 与 壮 板 便 度 分 别 呈 负 相 关 和 正 相 关 ， 表 明 随 
着 硬度 的 提高 ， 主 要 破坏 形式 由 弯曲 一 拉 伸 转变 为 剪 切 。400 HBW 由 于 发 生 骨 落 ， 残 
余 厚 度 较 低 。 各 尊 板 破坏 形式 均 与 临界 击 穿 相 似 ， 均 存在 环形 破坏 面 : 300 HBW 面 厚 
度 最 大 ， 背 表面 存在 损伤 但 未 完全 失效 ;400 HBW 发 生 背 面 月 落 ， 以 至 于 残余 靶 厚 反 
常 偏 低 : 500 和 600 HBW 靶 损 伤 形式 相同 ， 均 包括 环 状 剪 切 和 背面 弯曲 两 部 分 。 


表 3.$ 临界 阻挡 工 况 下 弹 靶 形 瑶 


300 HBW 400 HBW | 3S00 HBW | 600 HBVW 
| | | 员 必 了 
损伤 
人 
损伤 (复原 ) 
1 _ 
表 3.6 背面 最 大 变形 量 
AHBW Tsymm BFDMmm Fnmm 
30 9 571 471 
400 个 5. 55 3. 86 
500 6.5 5. 17 5 总 
600 6 4. 94 5. 10 
3.4 半 无 限 靶 


为 明确 半 流 体 侵 彻 性 能 ， 进 行 半 无 限 半 仿 真 。 取 靶 厚 25 mm 以 完全 排除 背面 变形 、 
应 力 波 反射 等 背面 效应 的 影响 。 由 图 3.3 可 见 ，300 和 400 HBW 靶 均 受到 显著 侵蚀 ， 弹 
体 碎 片 形 成 明显 管状 残余 体 ， 划 板 形成 翻 唇 ; 500 靶 翻 展 在 终 时 刻 已 断裂 ， 向 后 随 弹 体 
碎片 一 同 飞 散 ， 形 成 靶 前 碎片 云 ，600 HBW 靶 仅 受 轻微 侵蚀 。 如 图 3.5 所 示 ， 所 有 条 均 
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发 生 球 面 和 轴线 两 部 分 失效 : 300、400 和 5$00 HBW 间 的 球面 失效 区 域 均 深入 靶 体 内 部 ， 
在 其 外 部 存在 无 失效 区 ， 即 图 3.3 与 3.4 中 的 翻 唇 /大 碎片 ; 600 HBW 由 于 硬度 高 且 塑 性 
低 ， 球 面 失效 厚 而 浅 ， 不 存在 外 部 无 失效 区 ， 因 此 未 能 形成 翻 展 。 


(a) 300 HBW (b) 400 HB (oO500HBW (dj 600HBW 
图 3.3 终止 时 刻 弹 靶 形 貌 


| | 办 | 六 几 放 ee 
(a) 1.75 hs (b) 2.24 hs (c) 2.72 hs (d) 3.12 ns 
图 3.4 500 HBW 刘 翻 展 形 成、 断裂 与 靶 前 碎片 飞散 过 程 
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(a) 300 HBW (b) 400 HBW (o) 500HBW  ， (d) 600HBW 


图 3.5 终止 时 刻 有 体 损伤 


由 头 尾 速度 可 分 为 三 个 阶段 ， 以 下 以 侵 彻 深度 最 大 的 300 HBW 为 例 : 

1) 初始 瞬 态 阶段 : 约 0~0.5 mm。 此 阶段 弹 靶 开始 接触 。 由 于 弹 体 强度 远 低 于 靶 板 ， 
接触 面 上 的 弹 体 迅 速 变形 、 扩 张 ， 头 部 速度 迅速 下 降 

2) 准 定常 阶段 ， 对 应 侵 彻 深度 约 0.5~4.2 mm。 此 阶段 类 似 于 半 流 体 长 杆 侵 彻 的 主 
人 尾 速度 及 两 者 差 值 近似 保持 恒定 ， 其 中 头 部 速度 约 在 100 
m/s 左右 ， 尾 部 速度 接近 于 其 初始 速度 990 m/s。 根 据 实 际 发 生 侵 彻 的 部 件 ， 可 分 为 两 个 
小 阶段 : 


2.1 ， 被 甲 侵 蚀 阶段 ; 0.5~1.3 mm。 此 阶段 弹 体 强度 较 高 的 被 甲 部 分 与 半 体 接 
触 1.3 mm 时 尖端 被 甲 完 全 破坏 ， 头 部 速度 出 现 谷 值 ， 头 尾 速度 差 出 现 峰 值 ， 

2.2 ，” 铅 芯 侵蚀 阶段 ;1.3~4.2 mm。 此 阶段 为 铅 芯 弹 侵 彻 半 无 限 靶 的 最 主要 阶 
段 ， 铅 芯 在 此 阶段 与 间 板 在 轴 向 直接 接触 ， 被 甲 向 外 翻 折 ， 仅 对 铅 芯 有 微弱 的 径 向 
约束 作用 。 

3) 停止 阶段 : 4.2~5.5 mm， 相 当 于 半 流 体 长 杆 侵 彻 的 侵 彻 末期 。 此 阶段 弹 体 不 再 满 
足 流体 头 部 一 刚体 主体 假设 。。 可 分 为 三 个 小 阶段 。 

3.1 ”4.2~S.0 mm: 头 部 速度 仍 保持 在 约 100 my/s 不 变 ， 尾 部 速度 开始 下 降 。 
同时 靶 体 扩张 减缓 ， 至 5.0 mm 时 弹 体 直径 达到 最 大 值 ， 尾 部 减速 度 逐 渐 升 高 

3.2 ”5.0~5.5 mm: 尾部 快 减速 ， 阻 力 下 降 。 至 5.5 mm 时 头 部 速度 归 零 ， 刘 
体 不 再 受 破坏 。 

3.3 5.5~ 终 止 : 此 阶段 为 靶 体 的 弹性 回 弹 ， 弹 坑 深 度 略 微 减 小 ， 侵 彻 深度 随 
之 振荡 。 此 阶段 相对 于 整个 侵 彻 过 程 而 言 占 比 极 小 。 
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一 一 300HBW Head 
一 一 300HBW Tail 
一 一 400HBW Head 
一 一 400HBW Tail 
一 一 500HBW Head 
0 一 一 500HBW Tail 
一 一 600HBW Head 
一 一 600HBW Tail 


1000 5 


Velocity Difference (my s) 


Depth of Penetration (mm) 


(a) 头 尾 速度 


) 
六 
之 
| 


Resistance (kN 
信 
弓 
| 


Depth of Penetration (mm) 


(b) 头 尾 速度 差 


Depth of Penetration (mm) 


(c) 阻力 
图 3.4， MI193 侵 彻 各 硬度 钢 间 过 程 
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(al) DoP=0.3 mm (b) DoP=1.3 mm (c) DoP=4.2 mm (d) DoP=5S.0 mm 
图 3.5 300 HBW 弹 靶 变形 与 破坏 


尽管 各 钢 靶 均 可 观察 到 半 流 体 侵 彻 行为 ， 但 由 于 铬 必 弹 长 度 较 短 ， 且 不 由 单一 材料 
组 成 ， 其 侵 彻 过 程 中 速度 和 阻力 的 变化 与 典型 长 杆 弹 有 较 大 差别 。 由 图 3.4 可 见 ， 头 部 
速度 存在 较 大 幅度 的 波动 ， 整 体 上 与 靶 体 便 度 呈 负 相关 ， 同 时 头 尾 速度 差 与 间 体 硬度 正 
相关 ， 表 明 提 高 轩 板 硬度 能 更 有 效 地 损耗 弹 体 ; 阻力 在 整个 过 程 中 不 存在 稳定 阶段 ， 而 
是 始终 随 弹 体 的 钝 化 和 扩张 而 增 大 。 

由 图 3.5 可 见 ， 侵 彻 深度 与 最 大 阻力 均 与 靶 体 便 度 呈 近 似 线性 关系 。 与 3.1 节 对 最 
大 应 力 的 分 析 相 同 ，400 和 500 HBW 撮 的 抗 弹 能 力 分 别 较 拟 合 线 俩 低 和 侦 高 ， 表 明寺 
体 的 材料 模型 与 硬度 存在 误差 。 
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四 
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360 


Maximum Resistance Force (kN) 
CD 
己 
己 


二 全 仆仆 全 全 攻 生生 生生 全 用 全 生生 本 枚 站 全 全 交 人 全 向 全 二 
四 Simulation 
Fittin 

260 于 本 人 生生 生生 全 作 人 生 汪 生生 全 下 人生 和 仆 生生 人 且 全 和 汪 革 全 全 5 
240 T 1 
300 400 500 600 

Hardness (HBW) 

三 
(b) “最 大 阻力 


图 3.6 ”990 m/s M193 侵 彻 300~600 HBW 钢 靶 数值 结果 


3.5 ”等 效 系数 

分 别 以 对 半 无 限 靶 侵 彻 深度 CK7) 、 对 有 限 厚 靶 的 临界 穿 透 厚度 CKP) 衡量 相对 于 
300 HBW 靶 的 抗 弹 性 能 ， 并 最 大 应 力 (Ke) 辅助 比较 。Borvik 等 对 尖 头 刚体 弹 的 实验 
结果 表明 ， 苯 板 的 抗 弹性 能 仅 与 有 关 ， 并 提出 装甲 钢 的 强度 远 比 塑性 重要 01。 然 而 ， 由 
表 3.7 和 网 3.7 可 见 ， 三 种 等 效 方法 所 得 结果 有 很 大 偏差 : Kz 与 Ke 均 呈 现 与 硬度 的 良好 
相关 性 ， 表 明 对 于 仅 少 量 失效 的 半 无 限 靶 抗 弹性 能 与 壮 板 硬 度 /强度 高 度 相 关 ; KR 尽管 
同样 与 硬度 正 相 关 ， 但 400 HBW 明显 较 其 余 偏 低 ， 而 300、500 和 600 HBW 呈 与 另 两 
等 效 值 相似 的 线性 关系 。Kz 与 Kz 的 偏差 主要 来 源 于 靶 板 破坏 形式 的 差异 ， 对 于 临界 击 
穿 的 有 限 厚 靶 ， 其 塑性 决定 了 抗 拉 伸 和 剪 切 的 能 力 ， 其 中 对 于 强度 较 低 、 弯 曲 幅 度 较 大 
的 300 和 400 HBW 以 前 者 为 主 ， 对 于 强度 较 高 、 弯 曲 幅 度 较 小 的 5300 和 600 HBW 以 后 
者 为 主 。400 HBW 单 轴 拉 伸 失 效应 变 1.38, 双 轴 拉 伸 失效 应 变 0.40; 作为 对 比 , 300 HBW 
单 轴 拉 伸 失 效应 变 3.3$， 双 轴 拉 伸 失 效应 变 2.93。400 HBW 极 低 的 双 轴 拉 伸 失效 应 变 
导致 其 更 易 破坏 ， 因 此 需要 较 预 期 更 大 的 厚度 以 抵消 此 缺陷 的 负面 影响 。 


表 3.7 三 种 等 效 方法 比较 
PHBW Ar RE 的 
300 1 1 1 
400 “12 128 1.26 
500 1.38 2.07 1.64 
600 15 2.60 1.89 
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一 汪 一 Strength 
一 -一 thickness 


Effectiveness 
品 
坊 
CI 


300 400 500 600 
硬度 (HBW) 


图 3.7 ”三 种 等 效 方法 比较 


3. 6 本 章 小 结 

本 章 基于 2 节 优化 的 SPH-FEM 耦合 模型 ， 比 较 了 等 厚度 击 穿 、 弹 道 极限 和 半 无 限 
靶 三 种 靶 体 厚度 情况 下 铅 芯 弹 的 侵蚀 行为 。 结 果 表 明 : 铅 芯 弹 对 靶 板 的 破坏 包含 半 流体 
侵 彻 和 前 切 冲 塞 两 部 分 ， 其 中 中 低 硬度 肢 以 前 者 为 主 ， 高 硬度 靶 以 后 者 为 主 。 但 相 较 于 
生 型 半 流体 侵 彻 ， 铅 芯 弹 冲击 存在 非 均 质 弹 体 、 背 面 效应 等 非常 规 因素 ， 因 此 侵 彻 过 程 
中 的 速度 、 阻 力 等 量 的 变化 存在 差别 ， 而 相 较 于 典型 冲 塞 ， 铅 芯 弹 冲击 存在 弯曲 变形 、 
拉 伸 断 裂 等 异常 现象 ， 以 致 丢 板 的 毁伤 形 貌 与 理想 情况 不 同 。 
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4 理论 建 模 


综合 考虑 低 强 度 靶 的 半 流 体 侵 彻 和 高 强度 半 的 名 塞 两 种 破坏 形式 ， 本 文采 用 半 流 体 
侵 彻 体现 弹 靶 相互 侵蚀 的 过 程 ， 以 剪 切 诈 塞 判 定 击 穿 。 由 于 低 便 度 刘 的 破坏 除 剪 切 种 塞 
外 也 包括 拉 伸 ， 而 根据 Johnson-Cook 模型 ， 靶 板材 料 的 抗 拉 拨 性 显著 弱 于 抗 剪 ， 故 此 模 
型 求 得 的 对 低 硬 度 靶 的 极限 穿 透 厚度 会 较 实 际 值 往 低 。 此 外 ， 铅 忌 步 枪弹 几何 外 形 与 理 
想 钝 头 长 杆 有 较 大 区 别 ， 且 密度 在 轴 癌 和 径 向 上 分 布 缘 不 均匀 ， 不 能 使 用 简单 的 表达 式 
精确 描述 ， 因 此 本 模型 仅 用 于 定性 解释 和 分 析 侵 彻 过 程 。 

4.1 半 流 体 侵 彻 

Alekseevskii-Tate 半 流 体 侵 彻 模 型 提出 于 1967 年 55154, 首次 将 固体 的 高 速 侵 彻 与 流 
体 动力 学 相 结合 ， 是 长 杆 高 速 侵 彻 最 为 经 典 的 理论 模型 ， 也 是 后 世 绝 大 多 数理 论 模 型 的 
基础 。 由 于 A-T 模型 仅 体现 一 维 效应 ， 且 对 侵 彻 过 程 的 描述 过 于 平均 化 ， 其 计算 精度 较 
低 。 国 内 外 学 者 对 A-T 模型 进行 了 了 诸多 改进 ， 其 中 最 具有 代表 性 的 是 Rosenberg 和 
Walker-Anderson 模型 。Rosenberg 等 注意 到 侵 彻 过 程 中 弹 靶 界面 压力 从 弹 体 头 部 中 心 辐 
边缘 的 衰减 61。 为 体现 此 二 维 效应 ， 在 A-T 模型 的 基础 上 引入 刚性 杆 体 和 钝 化 头 部 / 轩 
体 等 效 截 面积 的 概念 。Rosenberg 模型 的 基本 形式 为 : 


1 1 
了 Po 一 六 十 也 -3GDmm 十 及 ) 


履 _ 

的 Po (4.D 
和 

帮 

人 

帮 


式 中 和 产 分 别 为 弹 体 和 苯 体 密度 ; ww 记分 别 为 弹 体 头 部 和 尾部 速度 ;了 乃 为 弹 体 强度 ， 
民 为 并 体 阻力 ;8$: 为 弹 、 划 等 效 截 面积 之 比 ; ! 为 弹 体 剩 余 长 度 ，P 为 侵 彻 深度 。 式 中 
厂 、R、S 需 由 经 验 估算 或 从 仿真 结果 拟 合 。$8 参考 Rosenberg 等 的 理论 分 析 575 和 王 儿 等 
5 的 实验 结果 ， 对 于 一 般 长 杆 弹 约 为 2; 在 A-T 及 其 它 不 考虑 二 维 效应 的 模型 中 为 1。 
四 个 方程 中 最 主要 的 为 首 个 方程 ， 通 党 也 被 称 为 平衡 方程 ， 用 于 计算 w; 后 三 个 方程 分 
别 计算 弹 体 减 速 、 弹 体 侵 蚀 和 间 体 侵蚀 。 由 于 铅 芯 和 被 甲 的 强度 均 远 低 于 四 种 装甲 钢 ， 
对 的 取 值 不 会 对 模型 求解 造成 显赫 影响 ， 初 设 为 30 MPa。 因 此 ， 模 型 中 的 未 知 量 为 等 
效 截面 积 比 % 和 郭 体 阻力 R。 
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及 可 通过 界面 击溃 现象 求解 ， 即 由 v=0 时 的 平衡 方程 一 一 式 (4.2) 一 一 求解 3$3 和 及。 
3.3 和 3.4 节 分 析 表 明 ,，300~500 HBW 和 600 HBW 部 的 损伤 形式 存在 明显 差异 ， 其 中 后 
者 在 临界 击 穿 时 塞 块 几乎 无 弯曲 ， 在 半 无 限 靶 时 不 产生 翻 恨 。 因 此 ， 当 靶 板 硬度 在 
500~600 HBW 之 间 的 某 一 数值 时 发 生 界 面 击溃 ， 即 侵蚀 速度 的 下 限 高 于 弹 体 初 速 990 
ID/S 。 


1 
5 局 记 = 及 (4.2) 


Rosenberg 模型 中 尺 由 式 (4.3) 计 算 。 式 中 c 为 靶 体 等 效 强度 ， 对 于 500 HBW 坦 为 
1637 MPa， 对 于 600 HBW 靶 为 1799 MPa; 玖 为 靶 体 弹性 模 量 ， 对 于 钢 为 206 GPa ， 
为 泊 松 比 , 对 于 钢 为 0.3。 由 以 上 参数 计算 , 则 500 和 600 HBW 的 玉 分 别 为 4772 和 $146 
MPa。 考 虑 到 500 HBW 靶 受 破坏 较 强 ， 且 600 HBW 靶 表 面 仍 有 少量 失效 ， 初 取 临 界 界 
面 击溃 时 RF=5000 MPa， 则 S=1.11， 显 著 小 于 钨 合金 长 杆 弹 侵 彻 钢 靶 的 等 效 截 面积 ， 表 
明 由 于 铬 芯 弹 速度 慢 、 强 度 极 低 ， 其 径 向 流动 、 侵 蚀 能 力 很 弱 。 

Or V3 已 
考虑 到 步枪 弹 的 几何 外 形 与 理想 的 圆柱 长 杆 有 较 大 区 别 ， 为 简化 计算 ， 将 原 弧 锥 一 

圆柱 的 铜 一 铅 结构 近似 为 全 铅 圆锥 一 圆柱 组 合体 。 对 于 M193 有 : 

0.761 ”0<1<7.5 

| 5$.7 7.5<17<19.3 

以 上 式 计算 990 m/s M193 侵 彻 半 无 限 靶 过 程 ，300~500 HBW 驾 侵 彻 深度 分 别 为 
3.79、2.84、0.86 mm， 准 定常 段 平 均 头 部 速度 放 分 别 为 162.2、127.1、42.6 m/s。 与 实 
际 侵 彻 半 无 限 靶 过 程 相 比 ， 理 论 计 算 所 得 侵 彻 深度 明显 偏 低 ; 平均 六 由 于 仿真 结果 存在 
大 幅 波动 ， 难 以 定量 判断 ， 但 其 范围 基本 准确 。 偏 差 主 要 来 源 于 弹 体 结构 的 简化 和 
Rosenberg 模型 自身 缺陷 。 受 限于 长 杆 半 流 体 侵 彻 研究 的 局 限 性 ， 目 前 尚 无 考虑 杆 体 非 
理想 结构 的 理论 计算 模型 。 尽 管 本 文中 铅 芯 和 被 甲 的 强度 尽管 均 远 低 于 靶 钢 ， 但 两 者 问 
仍 存在 很 大 差距 ， 因 此 平均 杆 体 强 度 的 选取 依赖 于 经 验 ， 不 能 通过 理论 计算 直接 获得 。 
张 清 源 的 对 照 性 研究 表明 ，Rosenberg 模型 本 身 预测 侵 彻 深度 偏 低 691。 两 者 的 共同 人 往 差 
导致 半 流 体 侵 彻 过 程 无 法 准确 计算 ， 但 侵 彻 深度 来 看 能 定性 反映 刘 板 的 抗 弹 能 


4.2 勇 切 冲 塞 


对 于 驾 体 几乎 不 受 侵蚀 ， 以 及 穿 甲 末期 靶 板 背面 显著 凸 出 ， 不 再 可 被 视 为 半 无 限 刘 
的 情况 ， 采 用 Recht 动量 一 能 量 平 衡 理论 判定 是 否 发 生 冲 赛 穿 孔 。Recht 理论 的 基本 假 


有 R = ] (4.3) 


(4.4) 
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设 为 : 1) 弹 体 和 靶 板 材料 不 同 ， 因 而 声波 波 速 不 同 ，2) 弹 体 不 是 刚性 的 ; 3) 弹 体 和 
塞 块 的 径 向 尺寸 可 能 不 同 。 高 速 销 心 弹 冲 击 钢 部 与 以 上 假设 均 相 符 。Recht 认为 ， 剖 塞 
穿孔 过 程 实际 上 是 一 个 动能 转换 的 过 程 ， 弹 体 原 有 动能 Bo 转化 为 : 

1) 弹 体 剩 余 动能 和 塞 块 动能 ， 即 前 文 的 罕 出 体 动能 Br; 

2) 弹 孔 周围 靶 体 材料 剪 切 破坏 吸收 的 应 变 能 现 ; 

3) 弹 驾 伴 撞 过 程 中 形成 共同 速度 消耗 的 能 量 到 。 
由 此 得 到 能 量 守恒 方程 : 


巨 。= 玖 ,+ 丙 + 丽 / (4.5) 
由 于 弹 体质 量 损失 和 新 体 的 压缩 变形 和 侵蚀 的 主要 过 程 已 在 半 流 体 侵 彻 部 分 计算 ， 
且 剪 切 破 坏 耗 时 相对 较 短 ， 故 假设 冲 塞 前 后 弹 体 质量 不 变 ， 塞 块 质量 由 剩余 靶 板 厚度 和 
弹 体 直径 确定 ; 球 仅 包括 剪 切 失效 应 变 能 。 因 剖 塞 发 生 时 弹 体 尖端 已 被 完全 匆 粗 ， 故 近 
似 弹 体 为 均匀 圆柱 体 。 由 此 得 到 实际 冲 塞 过 程 的 能 量 守恒 方程 : 

] 

9 
式 中 : mp、 几 分 别 为 弹 体 、 塞 块 质量 ， 可 由 半 流 体 侵 彻 模 型 粗略 确定 ，D, 为 塞 块 直径 ， 
参考 仿真 与 半 流 体 理 论 分 析 结 果 粗 略 计 算 时 取 弹 体 扩 张 后 直径 , 即 1.5 倍 弹 体 初始 直径 ; 
70、 六 为 弹 体 初 始 、 剩 余 速 度 ， 其 中 如 由 半 流 体 侵 彻 模 型 确定 ， 对 于 杆 体 〈 铬 芯 ) 强度 
远 低 于 靶 板 的 M193 可 取 初 速 作 为 近似 ; 志 为 靶 体 残余 厚度 ， 由 半 流 体 侵 彻 模 型 确定 ; 
六 为 剪 切 应 变 区 径 向 等 效 厚度 ， 与 靶 体 塑 忻 有 关 ， 需 由 仿真 结果 反 辐 推算 ;ws 为 剪 切 失 
效应 变 能 密度 ， 由 Johnson-Cook 材料 模型 确定 。 知 以 上 方程 有 对 六 的 正 实数 解 ， 则 判 
定 靶 板 被 击 穿 。 
以 990 m/s M193 临界 击 穿 6 mm Armox 600 靶 工 况 验 证 以 上 剪 切 冲 塞 模型 。 为 确保 
准确 性 ， 冲 塞 破 坏 初 步 以 更 精确 的 仿真 结果 计算 。 以 弹头 双 圆 驳 部 被 完全 侵 彻 作为 冲 塞 
判定 起 始 。 由 于 弹 体 尾部 速度 在 半 流 体 侵 彻 过 程 几乎 无 下 降 ,， 假设 弹 体 速度 w 在 冲 塞 开 
台 时 仍 为 990 my/s。 其 余 参 数 分 别 为 : 塞 块 与 弹 体 的 整体 速度 为 186 m/s, mp =2.4 g, 7 天 2.2 
g，DF8.6 mm， 太 4.8 mm， 解 得 r=1.0 mm。 尽 管 Recht 模型 理论 上 能 较为 准确 地 判定 
剪 切 冲 塞 ， 但 弹 体质 量 m 的 计算 依赖 于 较 精 确 的 几何 模型 。 知 选用 半 流 体 侵 彻 模 型 中 
简化 的 全 铬 柱 体 ， 则 上 述 工 况 中 的 六 将 增 大 至 2.6 mm， 较 实际 情况 显著 偏 高 。 


| 3 1Ip10 ， 
切 三 二 (1 十 1 )77 十 元 DD1O + 一 一 一 (4.0) 
站 。 ”” 27 十 711, 
忆 
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4.3 ”本章 小 结 

本 章 基 于 第 3 节 仿真 结果 ， 建 立 了 简化 的 Rosenberg 半 流 体 侵 彻 一 Recht 剪 切 冲 塞 
模型 。 受 限于 铅 芯 一 铜 被 弹 复 杂 的 结构 和 理论 模型 自身 的 简化 与 不 准确 性 ， 尽 管 近似 理 
论 模 型 可 定性 解释 其 穿 甲 原 理 ， 但 精确 判断 击 穿 仍 依赖 于 仿真 。 


和 
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结 论 


本 文通 过 FEM、EFG、SPH 三 种 方法 研究 了 高 速 铅 忌 步枪 弹 侵 彻 不 同 强度 钢 刘 的 行 
为 ， 探 索 了 初始 光滑 长 度 参 数 、 粒 子 近似 方法 和 失效 粒子 控制 方法 对 SPH-FEM 耦合 方 
法 仿真 精度 的 影响 ,基于 半 流 体 侵 彻 与 剪 切 冲 塞 模型 总 结 了 高 速 铅 芯 步枪 弹 侵 彻 钢 靶 的 
规律 ， 主 要 结论 如 下 : 

《1) 以 SPH 模拟 发 生 大 变形 和 失效 的 靶 板 核心 区 与 弹头 , 以 FEM 法 模拟 仅 小 变形 、 
受 固 定 约束 的 靶 板 外 周 区 的 SPH-FEM 耦合 法 能 最 准确 地 模拟 高 速 铅 芯 步枪 弹 侵 彻 钢 靶 
的 过 程 。FEM 和 EFG 法 受 限 于 网 格 ， 必 须 删 除 失效 单元 ， 因 此 仿真 过 程 中 质量 不 守恒 ， 
导致 侵 彻 阻力 偏 低 、 侵 彻 深度 偏 高 ， 且 不 能 准确 模拟 冲 塞 穿孔 过 程 。 

2) 初始 光 涓 长 度 参数 CSLH 为 1.4， 粒 子 近似 方法 为 重 整 化 近似 FORM 1， 人 失效 
粒子 控制 方法 为 部 分 失效 、 不 保留 体积 响应 IEROD 1 的 参数 组 合 能 最 有 效 地 抑制 SPH 
方法 拉 伸 不 稳定 的 缺陷 ， 且 具有 较 好 的 计算 稳定 性 和 效率 。 

《3) 不 同 强度 钢 间 受 高 速 铅 心 弹 侵 彻 的 过 程 基 本 相同 ， 主 要 分 为 半 流 体 侵 彻 、 背 表 
面 损伤 、 剪 切 冲 塞 三 部 分 。 侵 彻 初 期 弹头 速度 较 高 ， 与 靶 板 相互 侵蚀 ; 划 体 阻力 约 等 于 
其 抗 压 强度 ， 靶 体 受 力 面 积 约 为 弹 体 头 部 截面 积 的 两 倍 。 侵 彻 过 程 中 ， 划 板 背 表 面 受 整 
体 变形 影响 ， 主 要 承受 拉 应 力 ， 形 成 明显 隆起 并 发 生 拉 伸 损 伤 。 侵 彻 末 期 弹头 因 阻 力 而 
径 向 扩张 , 靶 板 受 侵蚀 而 厚度 减 小 , 形成 直径 与 弹 体 接近 的 柱 形 剪 切 破 坏 面 并 冲 出 塞 块 。 
中 低 强度 部 的 穿孔 以 半 流 体 侵 彻 为 主 ， 侵 彻 末 期 除 剪 切 冲 塞 以 外 伴 有 背 表面 拉 伸 破坏 ; 
高 强度 靶 板 的 穿孔 以 剪 切 冲 塞 为 主 ， 半 流体 侵 彻 识 度 相 对 较 低 且 与 抗 压 韧性 正 相 关 ， 剪 
切 破坏 面 较 明 显 ， 背 表面 仅 微 小 损伤 而 无 破坏 。 

《4) Rosenberg 半 流 体 侵 彻 一 Recht 剪 切 冲 塞 模型 可 定性 解释 铅 芯 弹 的 穿 甲 原理 ， 
但 受 限 于 射 弹 结构 的 复杂 性 和 模型 自身 的 不 准确 性 ， 精 确 、 定 量 判断 击 穿 与 否 仍 依赖 于 
仿真 。 
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附录 A kk 文件 主要 内 容 


本 节 为 经 第 2 节 优化 后 的 990 m/s M193 侵 彻 6 mm Armox 600 靶 K 文件 ,不 包含 过 
于 宛 长 的 节点 /粒子 (*NODE) 、 单 元 (*ELEMENT SPH 和 *ELEMENT _ SOLID ) 和 节 
点 集 (*SET NODE LIST) 关键 字 


*KEYWORD 
*CONTROL BULK VISCOSITY 
$# dl1 q2 type btype tstype 
1.5 0.06 -2 0 0 
*CONTROL CONTACT 
$ slsfac Twnpal islchk Shlthk penopt thkchg ofrien enmass 
0.3 0.0 1 1 1 0 2 0 
$ usrstr usrfric nsbcs interm Xpene Ssthk ecdt tiedprj 
0 0 0 0 0.0 0 0 0 
$ sfric dfric edc vfc 也 th_ sf pen _Sf unused 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 
$ jignore frceng Skiprwg outseg Spotstp ”Spotdel spothin unusedl 
尼 1 0 1 0 1 0.5 
$ isym nserod TwWgaps rwgdth TWksf icov SWIadf itho 储 
0 0 1 0.0 0.0 0 0.0 0 
$ shledg pstiff ithcnt tdcnof ftall unused2 
1 0 0 0 1 0.0 
*CONTROL HOURGLASS 
$ ihq qh unused1l 
9 0.05 
*CONTROL SPH 
$ ncbs boxid dt idim memory form Start ImaxavV 
1 01.00000E20 3 150 1 0.01.00000E18 
$ cont deriv ini isShow ierod icont lavVis isymp 
0 0 0 1 1 0 0 100 
$ ithk istab dl unused 


1 
*CONTROL TERMINATION 


$ endtim endcyc dtmin endeng endmas unused 
.00000E-5 0 0.001 0.0 100000.0 0 

*CONTROL TIMESTEP 

$ dtnit tssfac isdo tslimt dt2mas lctm erode mslst 


0.0 0.0 0 0.0 0.0 0 1 0 


二 起 必 二 
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$ dt2msf dt2mslc imscl unused 
0.0 0 0 0.0 
*DATABASE GLSTAT 
3# dt binary lcur ioopt 
6.00000E-8 3 0 1 
*DATABASE MATSUM 
8# dt binary lcur ioopt 
6.00000E-8 3 0 1 
*DATABASE BINARY D3PLOT 
4$# dt lcdt beam npltc psetid 
.00000E-7 0 0 0 0 
$4# ioopt 
0 
*DATABASE EXTENT _ BINARY 
4# neiph neips maxint strflg Sigflg epsflg rltflg engflg 
6 0 0 0 1 1 1 
4$4# cmpflg ieverp beamip dcomp Shge StSSZ n3thdt lalemat 
0 0 0 1 1 1 过 
4$4# nintsld ”pkp_sen Sclp hydro msscl therm intout nodout 
0 0 0.0 0 2 0 
*BOUNDARY SPC_SET 
$# nsid cid dofx dofy dofz dofrx dofry dofrz 
1 0 1 1 1 1 1 
*CONTACT TIED NODES TO_SURFACE 
$# SSid msid SSstyp mstyp Sboxid mboxid SPDT ImpT 
12 1 3 3 0 0 0 
$# 全 | dc VC vdc penchk bt dt 
0.0 0.0 0.0 0.0 10.0 0 0.01.00000E20 
$# sfs sfm SSt mst Sfst sfmt fsf VSf 
1.0 1.0 0.0 0.0 1.0 1.0 1.0 1.0 
$ Soft Softscl lcidab maxpar Sbopt depth bsort frcfrq 
2 0 0 0 3 5 0 
$ “penmax 人 khopt Shlthk Snlog isym i2d3d Sldthk Slds 寻 
0 0 0 0 0 0 0 
*PART 
3# title 
Componentl'out 
4$# pid Secid mid eosid hgid grav adpopt tmid 
1 1 1 1 0 0 0 
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*SECTION_SOLID 


aet 
0 
5 e PT dtf 
79230.0 206000.0 0.3 0.0 
卫 C Im tm 
0.173 ”0.00877 0.712 1537.85 
Spall 直 dl d2 
2.0 0.0 -0.4 1.5 
erod efmin numint 
01.00000E-6 0.0 
S1 S2 S3 gamao 
1.49 0.0 0.0 1.93 
mld eoSsid hgid grav 
1 1 0 0 
hmin hmax Sphini death 
0.2 3.0 0.01.00000E20 
mld eoSsid hgid grav 
3 3 0 0 
5 e PT dtf 
43890.0 ”11535000.0 0.31 0.0 
了 C Im tm 
0.42 0.009 1.68 9153.85 
Spall 直 dl d2 
2.0 0.0 0.0 2.604 
erod efmin numint 
01.00000E-6 0.0 


$# secid elform 
1 ] 
*MAT _JOHNSON_COOK 
$# mid TO 
17.8S000E-9 
$# a b 
1S80.0 9538.0 
$# cp pc 
4.320000E8 0.0 
8# d5 c2p 
0.0 0.0 
*EOS_GRUNEISEN 
$4# eeosid C 
1 ， 4370000 
$# V0 
0.0 
*PART 
3# 
in1l 
$# pid Secid 
12 9 
*SECTION SPH 
$# secid cslh 
9 1.4 
*PART 
3# 
jacket 
$# pid Secid 
13 
*MAT JOHNSON_COOK 
$# mid TO 
38.80000E-9 
$# a b 
111.69 04.09 
$# cp pc 
3.760000E8 0.0 
3# d5 c2p 
0.0 0.0 
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VPp rateop 
0.0 0.0 
位 epso 
22.0 1.0 
d3 d4 
-0.5 0.0 
a e0 
0.0 0.0 
title 
adpopt tmid 
0 
Start ”sphkern 
0.0 1 
title 
adpopt tmid 
0 
VPp Tateop 
0.0 0.0 
位 epso 
22.03.00000E-4 
d3 d4 
-0.02 0.28 


*EOS_GRUNEISEN 
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4# eeosid C S1 S2 S3 gamao 
3 3940000 1.49 0.0 0.0 1.97 
3# v0 
0.0 
*PART 
3# 
core 
$# pid Secid mid eosid hgid grav 
14 3 4 4 0 0 
*MAT JOHNSON_COOK 
$# mid TO g e PT dtf 
41.13400E-8 7000.0 ”20090.0 0.435 0.0 
3$# a b Dn C m tm 
1.0 03.0 0.103 0.41 1 320.85 
$ cp pc spall 让 d1 d2 
1.690000E8 0.0 2.0 0.0 2.0 1.77 
$ d5 c2p 
0.0 0.0 01.00000E-6 0.0 
*EOS_GRUNEISEN 
4# eeosid C S1 S2 S3 gamao 
4 2006000 1.429 0.0 0.0 2.74 
3# v0 
0.0 
*INITIAL _ VELOCITY _ GENERATION 
$zmnsid/pid Styp omega VX Vy VZ 
13 2 19314.8 0.0 0.0 990000.0 
$# XC ycC ZC DX ny DZ 
0.0 0.0 0.0 0.0 0.0 0.0 
$ Sid Styp omega VX Vy VZ 
14 2 19314.8 0.0 0.0 990000.0 
$ XC yc ZC DX ny DZ 
0.0 0.0 0.0 0.0 0.0 0.0 
*DEFINE COORDINATE SYSTEM 
$# cid XO y0 Z0 xl yl 
93 0.0 0.0 0.0 1.0 0.0 
3# xp yp Zp 
0.0 1.0 0.0 


a e0 
0.0 0.0 
title 
adpopt tmid 
0 
VPp rateop 
0.0 0.0 
妇 epso 
22.03.00000E-4 
d3 d4 
-3.4 0.0 
unused 
a e0 
0.0 0.0 
ivatn icid 
0 93 
phase irigid 
0 0 
ivatn icid 
0 93 
phase “unused2 
0 0 
Z1 cidl 
0.0 0 
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修改 记录 


一 、 毕 业 论文 〈 设 计 ) 内 容重 要 修改 记录 
指导 教师 : 

1) 扩充 理论 建 模 章 

2) 调整 附录 A 格式 
评阅 教师 : 

1) 修改 公式 缩 进 格式 

2) 调整 图 表格 式 
答辩 委员 会 : 
1) 扩充 摘要 
2) 修改 页 眉 格 式 
二 、 毕 业 论文 〈 设 计 ) 外 文 翻译 修改 记录 
1) 修改 图 表 标 题 格式 

2) 调整 标题 分 级 格式 
三 、 毕 业 论文 〈 设 计 ) 正式 检测 重复 比 
总 文字 复制 比 1.7%， 重 复 字 数 529 


记录 人 (签字 ) , 王 息 锐 
指导 教师 〈 答 字 ) ， 黄 - 疙 丰 
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